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We study out of equilibrium dynamics and aging for a particle diffusing in one dimensional 
environments, such as the random force Sinai model, as a toy model for low dimensional systems. 
We study fluctuations of two times {tw,t) quantities from the probability distribution Q{z,t,tw) of 
the relative displacement z = x{t) — x{t^) in the limit of large waiting time t„ — > oo using numerical 
and analytical techniques. We find three generic large time regimes: (i) a quasi-equilibrium regime 
(finite t — t — t^) where Q{z,t) satisfies a general FDT equation (ii) an asymptotic diffusion 
regime for large time separation where Q{z)dz ~ Q[L{t) / L{ti^)]dz / L{t) (iii) an intermediate "aging" 
regime for intermediate time separation {h{t)/h{tjn) finite), with Q{z,t,t') = f{z,h{t)/h(t')). In 
the unbiased Sinai model we find numerical evidence for regime (i) and (ii), and for (iii) with 
^ , Q{z,t,t') = Qo{z)f(h{t)/h{t')) and h{t) ~ Int. Since h(t) ~ L(t) in Sinai's model there is a 

O ' singularity in the diffusion regime to allow for regime (iii). A directed model, related to the biased 

Sinai model is solved and shows (ii) and (iii) with strong non self-averaging properties. Similarities 
and differences with mean field results are discussed. A general approach using scaling of next 
highest encountered barriers is proposed to predict aging properties, h{t) and f{x) in landscapes 
, with fast growing barriers. It accounts qualitatively for aging in Sinai's model. We also identify 

^ ' a mecanism for aging in low dimensional phase space corresponding to an almost degeneracy of 

barriers. We illustrate this mechanism by introducing a new exactly solvable model, with barriers 
. and wells, which shows clearly diffusion and aging regimes with a rich variety of functions h{t). 

B 

a 
o 

O , I. INTRODUCTION 

There is presently considerable interest in out of equilibrium dynamical processes. For systems without quenched 
I disorder these are important to understand phenomena such as coarsening and domain growth^. In systems such as 
~V . spin glasses'*"*, random fields, interfaces, glasses in vortex systems*^'*^'*^, which are dominated by quenched disorder 
and ultra slow relaxations, a detailed understanding of out of equilibrium dynamics becomes absolutely necessary to 
ly-j ,, make contact with numerical simulations and experiments. These usually involve studying relaxation dynamics from 
■ an initial configuration at i = (e.g uncorrelated) and asking about correlations in the systems between two later 
C — times t' (also called the waiting time) and t < t' . An interesting question to ask is what happens when both t and 
t' are taken to infinity. Since there are of course many ways to take t, t' to infinity and one wishes to classify the 
possible regimes. 

In out of equilibrium situations the usual properties of equilibrium dynamics do not hold. Such properties are 
the time translational invariance (TTI) , i.e the dependence on t — t' only of the correlation functions as well as the 
fiuctuation dissipation theorem (FDT) which relates linear response to time derivatives of correlations functions. A 
(— I first question is then to ask how to take t, t' to infinity and still recover an equilibrium regime. 
Q Glassy systems with quenched disorder were found to exhibit a variety of non equilibrium properties often generically 

O ' termed "aging". Loosely it means that the properties of the system are governed by the age of the system t^, i.e the 
^ I time after the quench*^'**'*^ For instance it is expected that correlation functions in these systems have dependences 
such has t/tw This type of dependence is also found in simpler out of equilibrium systems such as coarsening in 
spin systems without disorder'^ which also exhibit dependences of correlation functions of the form L{t) / L{tw) ■ This 
dependence originates from the growth of domains of size L{t) and diffusion of the domain walls. A stronger form 
of aging seems to be observed^'^, and was proposed for spin glasses where the linear response shows memory effects, 
e.g the remanent linear magnetization after applying a field during time decays very slowly over a time scale 
set only by Other puzzling phenomena such as memory under thermal cycling are observed^'^'^. There is at 
present no theory which would account fully for all these phenomena^, and understanding is only partial. Ideas and 
scaling arguments borrowed from domain growth and coarsening were also applied to disordered models ("droplet 
picture"®), but it is unclear whether they can account for all situations. More recently exact solutions were obtained*^ 
for the out of equilibrium dynamics of several mean field models^^'^^ , some of which where found to exhibit strong 
aging properties. It makes several non trivial predictions*^ for both the correlation function C{t,t') and the response 
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function R(t, t') of these mean-field models in the limit of large times t.t' . The asymptotic time regime in these models 
has been successfully resolved, under some physical assumptions, and turns out to be in direct correspondence with 
Parisi's static replica symmetry breaking solution (with some important differences). However the matching of the 
small time regime to the asymptotic one remains problematic. A general picture for aging dynamics in mean field 
was proposed^^. The resemblance with some of the features observed in spin glass experiments seems encouraging, 
though many points remain unclear^". There has been some attempts to classify the various aging behaviours, and to 
differentiate the domain growth coarsening type of "aging" from a stronger type analogous to what is found in mean 
field models^°. Despite these recent advances, there is however at present no detailed microscopic understanding of 
the aging phenomena. 

An elegant microscopic mechanism for aging was proposed some time ago by Feigelman and Vinokur (FV)^'^ in 
the context of diffusion in a one dimensional environment. Using a semi-quantitative analysis, they proposed that 
traps with wide distribution of waiting times (a diverging first moment < r >= oo ) would naturally lead to aging 
phenomena and waiting time dependence. The idea is that because < t >= oo, at time the system is typically 
in a trap of release time ~ and thus the diffusing particle sees potential barriers which effectively grows with tyj. 
The Feigelman- Vinokur trap model was later popularized by Bouchaud^^ and generalized, on the basis of previous 
work on wide distributions of waiting times^*~^°'^^'^^ to describe various situations, some inspired by the statics of 
mean field models^^. It is important to note, however, that in the one dimensional Sinai model with a bias originally 
studied by FV and in Ref.^^ the wide distributions can be shown to be dynamically generated and not artificially put 
by hand. Despite the phenomenological appeal of the FV mecanism for aging it is unclear how far it can be pushed^^ 
to describe all the physics of glasses and provide a non artificial, dynamically generated mecanism for aging. In fact 
rather different scenarios were proposed to understand more microscopically mean field dynamics^^ 

In this paper we study some simple one dimensional diffusion model, such as Sinai's model and investigate in some 
details their out of equilibrium dynamics. We will not attempt to propose any phenomenological model for aging 
in real glasses but simply study the extreme case of low dimensional phase space and identify the various possible 
large time regimes. These low dimensional models are dominated by activated dynamics over energy barriers and 
stand at the opposite end from mean field models. One can hope that they serve as toy models for diffusion in low 
dimensional space and give some insight into finite dimension. In particular, one can check whether some of the ideas 
introduced in mean-field theory carry through. Specifically, mean field generalizations of Sinai's model (i.e diffusion 
in a dimensional random potential with d oo) were solved^^ and can be directly compared. The advantage of this 
one dimensional model is that the correlation and response functions can in principle be computed numerically up 
to very large times. Note that a simulation of Sinai's model was performed recently in^**, but with a much too slow 
algorithm to reach significant times, as we will discuss. A first summary of the present study was contained in^^. 

It becomes clear in our study that here the full probability distribution should be studied, while in mean field theory 
it is enough to study the second moment. For a diffusing particle the interesting quantity is the distribution Q{z, t, t') 
of displacements z = x{t) — x{t') between times t and t' (the particle having started at i = at a random position) 
with a translational or configurational average of Q. The second moment B{t,t') =< >q was computed in mean 
field^^. The results are as follows for t ^ oo and t' ^ oo. There is: 

(i) a quasi equilibrium regime for finite t = t — t' where TTI {B{t,t') = B{t)) and FDT theorems hold. In that 
regime displacements arc bounded since B{t ^ oo) = 60 • 

(ii) for more separated times B{t,t') > bo keeps growing. This is the aging regime where B{t,t') remains a fixed 
number B{t,t') = B provided t and t' grow in some well defined way, i.e with h{t)/h{t') a fixed number (a function 
of B). In that very non trivial regime some new generalized FDT theorems hold, according to the general theory 
of^^. The function h{t) is not determined by the large time mean field solution, and at present must be determined 
numerically. The number B can eventually be chosen as large as wanted. For large B this regimes crosses over into a 
diffusion regime. Interestingly there is a singularity at the beginning of the aging regime with a non trivial exponent 
/3, i.e the function t') = B[h{t)/h{t')] is non analytic as a function of h{t)/h{t') — 1. A singularity at the beginning 
of the aging regime is indeed found in a variety of experimental glassy systems. The exponent /?, relevant for mode 
coupling theories of real glasses^^, was obtained analytically in^^. 

Note also that in mean field there is some mathematical correspondance between two times t, t' and replica pairs a, b. 
There, roughly speaking the quasi equilibrium regime (i) was found to correspond to the replica symmetric part of the 
static solution, while the aging regime (ii) corresponded to the RSB part of the static solution. The correspondence 
is in fact not perfect because dynamical quantities do not always coincide with their static counterpart^^'^^'^^. 

It is important to know what remains of the above dynamical mean field scenario in low dimensional models, such as 
the models studied here, and we will attempt to give some elements of answer. Obviously this scenario will be modified 
since we now have to deal with full distributions Q{z, t, t'). One would like e.g. to identify possible mechanisms which 
can reproduce aging with non trivial functions h{t) and exponents (3. 

We will use Sinai's model as a starting point. It is interesting because it is related to coarsening models with disorder. 
For instance it can model a single interface in a random field model, or motion of kinks in vortex lines, or dislocation 
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loops, in presence of point disorder. The dyn amics of S inai's model has been studied extensively^^'^'*'^''''^'''^^. The 
model was shown to exhibit ultra-slow diffusion < x{t)'^ > ^ (log^)**. The response to a driving force / was also shown 
to be quite anomalous, with sev eral phases . There is a threshold force such that for f < fc the velocity vanishes 
V = and diffusion is sublinear < x{t)'^ > ~ with a continuously variable exponent /i = f / fc- For f > fc one 
has V {f — fc) and similar transitions in higher order moments of the displacement take place at larger /. The 
physics of this anomalous response was also understood and shown to be related^^'^^ to the existence of barriers 
with an exponential distribution of energy heights, resulting in power-law distributions of trapping times and Levy 
distributions for first passage times. 

Let us illustrate some of the questions by some simple consideration. In the case of the case of an applied bias, 
even the first moment of the displacem ent already contains interesting information. It is known^**'^^'^^ that one has 
exactly, at large times < x{t) — x{0) > = Ct^ where C depends on the details of the model. The aging nature of this 
expression is clear since one can write: 



< x{t) - x{t') > = C{t^ - t't") ~ C^i^ (1) 



tu 



for t' = tyj, t = tyj + T. This ind icates that there is an aging regime in this model since one can take t ^ oo and 
— > oo keeping < x{t) — x{t') > a fixed number, provided: 



(2) 

In particular, in the limit where tu, goes to infinity first, there is no motion at all in any finite time interval! We also 
note that the above result can be rewritten as: 

hit) 



<x{t)-x{t')> = C\n-^^ (3) 

where h{t) = cxp(f^). This form (3) (with h{t) unspecified) is of the type found in large time solutions of mean 
field models as mentionned above. Similarly, subaging forms for h{t) = exp((lnf)") with a > 1, as observed in 
experiments^'^, could in principle originate from logarithmic diffusion process. Indeed one can write: 

ln'*f-ln"f'=ln^- (4) 

with the observed h{t). One would thus like to investigate these analogies further. 

We will thus first study Sinai's model in Section II (symmetric model) and in Section III (directed model). Then 
we will give in Section IV a general discussion of barrier mechanisms for aging which we believe allow to understand 
the results obtained on Sinai's model. We introduce a method to estimate aging functions by looking at the sequences 
of next highest barriers. Finally, in Section V we introduce and solve a model based on these considerations which 
exhibits non trivial againg and diffusion behaviour. 



II. SYMMETRIC SINAI MODEL 



A. the model 

In this Section we study the one dimensional Sinai modeP^'^^'^^'^^'^'''^^. In its continuous version it is described by 
the Langevin equation: 

(It 

= F{x{t))+n{t) (5) 



where < r]{t)r]{t') >= 2T5{t — t') is the thermal noise and F{x)F{x') = a6{x~x') is a quenched random force which 
is gaussian and uncorrelated. It may have an average F{x) = f. Writing F{x) = —dU{x)/dx it describes the thermal 
motion of a particle in a one dimensional gaussian random potential landscape U{x). The landscape {x,U{x)) can 
be seen itself as a trajectory of a random walker [x playing then the role of the time). In the general case / > 
the lanscape is tilted (the walker experiences a bias) but we will only consider here the case / = (symmetric Sinai 
model). Then the random potential has long range correlations: 

{U{x) - U{y))^ = a\x - y\ (6) 
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Thus barriers grow with scale which results in an anomalously slow dynamics < {x{t) — a;(0))^ >~ (Int)^. The N 
dimensional version of this model (i.e x = (xi, X2, , x^) ) was solved in the limit TV ^ oo in Ref.'^^. 

In this Section we study a space discretizcd version of Sinai's model which can be easily studied numerically. It is 
expected that this discretized version has the same large time physics as the continuous model (this was shown for 
some quantities in the weak disorder limit^^'^^. It is defined by a Fokker-Planck (FP) equation: 

= {HFp)nmPm{t) = W„,„+lP„+l (i) + W^„,„_lP„_l (t) - + Wn-l,n)Pn{t) (7) 

where Pn{t) is the probability that a particle is at site n at time t (with some initial condition at t = 0). The 
hopping rates W„+i,„, W„_i,„ are quenched random variables. It is convenient to parametrize these rates as in Ref.^^ 
in the form: 

Wn-l,n = Wn,n-l=e'^- (8) 

This describes effectively the Arrhenius diffusion of a particle; on a one dimensional lattice in a random potential 
Un = —'^J2k=o 'i^k- There is in effect a force 20„ on the link between site n — 1 and n. The equilibrium solution of 
equation 7 corresponding to zero link current Jn,n-i = Wn,n-iPn-i — Wn-i,nPn = is P^' = e~^". Here temperature 
is set to T = I. In the discrete Sinai model the variables 0„ are chosen gaussian independent from site to site, with 
4>n4>n' = cr/4J„rt' and the random potential thus follows a gaussian discrete random walk as a function of n, with 
correlations {Un — Un')"^ = cr|n — n'|. 

The important quantity to determine is the Green function P(n, f|no,to) {t > to) of the FP operator Hpp in a 
given environment. It is defined as the solution of (7) with initial condition P(n,to\no,to) = Snno- 

We have computed numerically P(n, t\no, to) using exact diagonalization of the corresponding Schrodinger operator 
on a finite size ring. This operator is a tridiagonal symmetric matrix which can be easily diagonalized for large size L. 
The method, as well as the expressions for the correlation and response functions, is detailed in the Appendix A. We 
use L+1 sites k = 0,L with both reflexive and periodic boundary conditions. We have used up to L = 250 sites and 
averages over 10^ and 5 10^ disorder configurations. Because Sinai's diffusion is so slow, we were able to study times 
up to 10^^ without spurious effects (edges, precision). We used cr = 2 in all simulations and checked the consistency 
of the results with several random number generators. 

From this exact expression one can average either one or a product of two of these Green functions to obtain single 
time or two times quantities, respectively. For finite times, translational averages (with e.g uniform initial measure 
in an infinite single environment) and averages over disorder should coincide. Note that this will not be correct in a 
very special case studied below of a periodic environment. We now describe the results and their interpretation using 
simple arguments. 



B. single time quantities 



Averaged single time quantities can be obtained from the disorder averaged Green function P(n, t\no, 0). There are 
some exactly known results for this single time quantities and we will thus start by comparing with these results (as 
a check of our simulations). 

In^* the averaged probability density at the origin P(no, t\no, 0) was computed exactly in the continuum limit for all 
times t. It describes the weak disorder universal behaviour (see^^). At large t the complete scaling form for the averaged 
diffusion front, i.e the distribution of the variable n/{logt)^ was obtained in^'' as P{n,t\0,0) ~ (logf)~^p[n(logt)~^]. 
The function p{x) is (up to a constant rescaling): 

K-) = -E^-K-^-^N) (9) 

fe=0 
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FIG. 1. Single time averaged diffusion front: at large time our simulations cannot be distinguished from Kesten's prediction. 
We denote x = x/xrms{i). L = 125 and 5 10^ configurations. 
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FIG. 2. Same as Fig. 1. in logarithmic scale 
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Wc have computed numerically P{n,t\0,0) and the x ^ (logt)^ law as a check of our program. In Fig. 1 and Fig. 
2 we have plotted the scaling function p[x] as it is determined numerically and the analytic result of Kesten (wc have 



fixed the scale by imposing equality of the second moment and we denote Xrms = < x^(t) > ) . One can see that 
the agreement is excellent. It improves considerably on an earlier determination by Nauenberg*^. The behaviour 
X ^ (log(t/to))^ is plotted in Fig. 3 and also agrees. We have not attempted to fit the amplitude to known results 
(we have checked that the order of magnitude is correct) since we were satisfied with the agreement in Fig. 1. 
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FIG. 3. Numerical data compared to Sinai's diffusion law x ~ Xrms(t) ^ log{t/to)^. L = 125 and 5 10^ configurations. 
We now turn to two time quantities. 



There are no analytical results available at present, to our knowledge for two time quantities. This make the 
numerical simulation all the more important. 

We have computed numerically the two time averaged Green function: P{n,t\n' ,t')P{n' ,t'\no,to = 0). It is a 
complicated object which is difficult to analyze. Thus we have started by computing, as in mean field models, the 
configurationally averaged mean squared displacement. 



The best way is to plot it as a function of t — t' for different waiting times tyj = t'. It is represented in Fig. 4 One 
sees that it clearly does not depend only on t — t' and that, the larger the waiting time, the slower it grows. The 
dynamics slows down considerably as ^ oo. This is very reminiscent of what happens in mean field models. One 
also sees that an effective plateau develops as t' becomes larger. In mean field it would be a true plateau"^^ which 
defines the Edwards Anderson order parameter bsA- Thus here an important question is whether there is a finite 
limit for B{tw + T,tw) as tw ^ t ^ 1 as in the mean field model. 



C. two time quantities 



B{t, t') = {{x{t) - x{t')Y) = Y^{n- n'fP{n, t\n', t')P{n', t'\nQ, to = 0) 



(10) 
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B(t, t-x ) 
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FIG. 4. mean squared separation B{t = tw + t, tw) as a function of r for several t ^ tw A limit curve B{t) can be seen, but 
B{t) keeps growing with r while in mean field it would go to a constant bsA = B{t = oo). L = 125 and 10^ configurations. 

The answer to this question is negative. This is because, as will be discussed below, in finite dimension contrary 

to mean field, while the distribution of x{t) — x(t') converges towards a limit, it second cumulant is infinite. In fact 
studying only the second cumulant B{tw + r, t^) is of little help in the present problem. One needs to study the full 
disorder averaged distribution of relative displacements: 



Q{z,t,t^)= j dx'P{x' + z,t\x',t^a)P{x',tyj\0,to = 0) (11) 

and not only its second cumulant. We have found that when looking at that distribution, several regimes can be 
identified. 



D. quasi-equilibrium regime 

The first question to ask is whether there is a quasi equilibrium regime when t and <^ are close together. It can be 
defined, for instance as t — tw + t with r fixed and finite (and tw oo). In mean field this would correspond to the 
FDT regime, where time translational invariance and the fiuctuation dissipation theorem are found to hold. 

One can thus ask whether the following limit distribution exists: 

Q{z,t) = lim lim Q{z,tu, + T,tw) (12) 

tyj — ^oo L — ^oo 

where L is the size of the system and whether in that limit, some equilibrium theorems hold. 

We are not able to answer to this important question rigorously in all generality and we encourage other workers 
to do so. We will however provide some elements which we hope will shed some light on the issues. 

One strategy is to first study the case of an infinite periodic medium but with a very large period Lq (i.e almost 
disordered). There some things can be shown when tw is very large. Of course this is cheating a bit since the idea 
would be to take Lq ^ oo eventually, and thus this is like interverting the limit — > oo and Lq ^ oo. It does give 
some insight though. Thus we will then check, using numerics and physical arguments, whether it can be extended 
to a non periodic case. 

1. a periodic model 

We are considering an infinite periodic medium. We assume that the potential U (x) is periodic which corresponds 
to the situation with no drift. One can generalize to the case with a bias by taking a periodic force (using the results 
of*6. 'We will first consider a single environment (no configurational averages). 
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We are interested in the distribution: 

r+oo 



/-l-oo 
dx'PLoix' + z,tyj +T|x'i^|a;oO) 
-oo 



(13) 



It can be written as a sum over periods: 

QLo{z,tyj+T,ty,) / dx'PLoix' + k + Z,tyj + T\x' + k,tyj)PLo{x' + k,t^\xoO) (14) 

k Jo 

= [ dx'PLoix' + Z,T\x',0)PLo{x',ty,\xoO) (15) 
Jo 

where we have defined the periodized Green function: 

PLoix',t'\xoO) = J2Pi^oix' + kL,t'\xoO) (16) 

k 

and wc have used the periodicity of the disorder PLoix,t\x' ,t') = PLoix + mL,t\x' + mL,t'). So, up to now it is 
exact for any configuration of the disorder. 

Since the periodized function P is in fact the Green function of the problem on a ring of size Lq (see e.g^^'^^), it 
converges towards the equihbrium Gibbs measure on the ring in the large tw limit: 

g-U(x')/T 

lim PLoix', t^\xoO) = J dx'e-u(-'yT (^^^ 



This implies that: 



QLoiz,T)^ [ ° dx'PLoix' + Z,t\x'0)-^^ 

Jo J " 



dye-^iy)/^ 



(18) 



This is particularly useful when r is fixed and much smaller than the time necessary to travel a period (i.e (Inr)^ <C 
Lq in Sinai's model). This formula can also be generalized in presence of a drift using the stationary distribution 
obtained in^^. Note that these arguments extend to any finite dimensional problem. 



S. non periodic case 



In the non periodic case, when Lq oo before tw, it is obvious that for a single configuration of disorder there is no 
limit to Qiz, +t, <^) as — !■ oo. Indeed numerically Sinai's diffusion consists of sudden jumps to deeper and deeper 
wells. Thus even though in each successive well there is prcsmnably a quasi-equilibrium regime (3(z,t^ + r, t^,) ~ 
Qweniz,T) it will depend on the details of each new well encountered. In some sense there is a distribution of such 
quasi equilibrium distributions (which is not unlike the image from replica symmetry breaking). In addition, there 
are times where the packet jumps, but they presumably become more and more ra re at la rge t^- Thus only the 
disorder averaged (or translational average) can be expected to converge to some limit Qiz,T), as these features will 
get smoothed out. We will assume that this convergence holds (based on our numerical evidence). 

One possible further assumption is that once averaged over disorder one has: 

Qiz,T) = Qeqiz,T)= lim QLoiz,T) (19) 

I/O— »oo 

where QLoiz,T) is the distribution in the periodized medium (using the same disorder distribution) given in (18). 
We have called Qeqiz,T) this distribution. 

A possible reasonable starting assumption is that they are equal in Q = Q^q in Sinai's model. We do have good 
numerical indication that this is indeed correct, as shown in Fig. 5. However, because differences, if they exist, could 
be subtle, small and hard to detect we emphasize that whether this is strictly correct should be checked further (there 
are indeed cases where it is wrong - see the solvable model of Section V). In any case, the above equation provides at 
least a basis for comparison. 
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Q(z,T^le3, tw) 
le-01 




FIG. 5. Convergence, for a fixed r, of Q{z,tn, +7", t„) (represented here for = l(f and = 4 10^^) towards Qeq{z,T) 
computed here from formula (18,19) (curve called "t^ = infinity" in the figure). The limit curve Q{z,t = oo) is shown by 
comparison (curve called "fit"). L = 200 and 5 10^ configurations. 



3. correspondence with the statics 



Let us investigate the large r limit of Q{z,t). One expect, if the above is assumption is correct, that: 



g-J/(x'+z)/Tg-!7(x')/T 



/"^o 

Q{z,t) lim Qlo{z,t)^ Yim Ql^{z)= \im / dx' '^ -um/t\2 (^0) 

io^oo i/o->oo Lo^qoJq [J dye '-'\y>i^y 

Thus the two time calculation becomes a two replica correlation function in a replica calculation. 

properties of the quasi static distribution 

For Sinai model one can argue that there is indeed a limit averaged quasi static distribution Q{z). Prom the above 
formula (20) and since there is typically one global minimum which dominates the period Lq in Sinai's landscape, it 
is clear that this distribution will consist of two parts: 

(i) a part localized near z = of finite extent, with a finite weight. 

(ii) a part which is strongly non self averaging with, for each environment, some localized peaks around some Zi. 
Averaging over disorder yields a smooth Q{z) with algebraic tail at large z: 

Qi^) ~ ^ (21) 

This tail comes from configurations of the disorder which look as in Fig. 6. The lowest well Umin in the sample of 

period Lq is represented and it happens that there is a secondary well, with a bottom at Umin + En with En = 0{T), 
at a finite distance zq of the first one. Then roughly the measure Q{z) will consist of two peaks localized around the 
two wells. 

Since in Sinai model the random potential can itself be seen as performing a random walk, the probability that such 
an environment occurs can be estimated from the probability of return to the origin of a random walk. This yields a 
probability ^ Zq^^'^ and thus the above algebraic tail of Q{z). More refined calculation, taking into account that the 
principal well is an absolute minimum can be performed and lead to estimates for A, but go beyond this paper. Let 
us point out that the above equilibrium Gibbs measure has been analyzed recently in Ref.^^ and they also obtained 
analytically the above algebraic tail (i.e formula (21)). This tail however can be explained from simple arguments. 

It would be interesting to check whether one can extend ideas from mean field and whether in a replica calculation 
of (20) the localized part near the origin would correspond to the replica symmetric part of the solution while the 
tails from the rare events to the part with broken replica symmetry. 
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So there is also an Edwards Anderson order parameter here but it is a distribution. Because of non self averaging 
effects one needs the full distribution Q{z). It would be nice to check more rigorously that Q{z) — Qeq(z) (and thus 
that there is no extra dynamical order parameter). 

Finally let us emphasize that the above arguments for the tail of the quasi equilibrium distribution are very similar 
to static droplet model arguments^. Here however they are made in a dynamical context. 



U 




FIG. 6. Two well (FDT-TTI) quasi equilibrium regime: rare configurations leading to sample fluctuations and algebraic 
decay of Q{z) 



The above arguments thus suggest that the moments of the relative displacements: 

lim < {x{t^ + t) - x{t^))^ > ~ / z"Q{z,t) 



(22) 



converge towards a finite limit when r — > oo only for n < 1/2. For n > 1/2 they grow unboundedly with r. These 

moments growing with t may not a priori be incompatible with being in a quasi equilibrium state (though it calls for 
further investigations). In mean field, this does not happen since even for unbounded diffusion problems'^® the quasi 
equilibrium regime is such that B{t, t') < hsA- 

We now estimate these moments (22) using a simple two well model taking into account the crossing of the barrier 
El, between the two wells. Let us write: 



Q{z,t) = 6{z)^{l + e— ""^) + 5{z - ^o)^(l - e— ""^) 



(23) 



1 /2 

and average over zq and the barrier height which we can take to scale as Ei,{zo) ~ azQ a being a positive random 
variable. This yields the moments: 



lim < {xity, + r) - a;(f,„))" > ~ ^ / dzz"^(l - e-^"'"''') 

tyj ^OO Z J Z ' 

Introducing the scaled variable y = z/ In^ r and using that the function (1 — e~'' ' 
finds, for n > 1/2. 



(24) 



) w BiXjo? - y) one 



lim < (a;(t^ + t) - a;(i^))" > ~ A(lnT) 



1n-\ 



(25) 
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FIG. 7. Numerical determination of the moments < |a|"(T) >q as a function of Inr. L = 200 and 5 10^ configurations. 

We have plotted these moments as determined numerically in Fig. 7. The logarithmic growth of moments n > 1/2 
is clearly demonstrated (we have checked that e.g. n = 1/4 saturates at large r). Though they seem to follow (25) 
qualitatively, a more quantitative agreement is probably difficult to reach numerically - since these effects entirely 
come from rare events ! 

The above considerations also predict that the distribution Q[z,t) should converge towards Q{z) with a z/\r^T 
scaling behaviour. This is consistent with our simulations as can be seen in Fig. 8 where we have shown the Qeq{z, r) 
for various r. It does indeed converge towards a limit curve Qeq{z) which was well fitted by its asymptotic behaviour 

^ ,,-3/2 

One can also check from this simple model (23) that as a consequence, the Edward Anderson "overlap" parameter 
Q{z = O.t) should also converge towards it equilibrium value as r ^ oo with l/(lnT)" corrections. This indeed 
happens as is shown in Fig. 9. 
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FIG. 9. Plot of Qeq{z = 0, t) — Qeq{z = 0,T = oo) against t (in Inr scale). L = 200 and 5 10^ configurations. 



4.. equilibrium theorems 

It is plausible that for large tu, equilibrium theorems like FDT hold for disorder averaged quantities. These provide 
relations between correlation functions and linear response functions. 

It is useful to generalize these theorems to the full probability distribution Q{z,t,t'). This is done in the Appendix 
A and B to which we refer for details. One defines: 



SQhiz,t,t') 



whore Qh{z,t,t') is the probability when an additional infinitesimal uniform field pulse h{t) = hS(t — t') is applied 
at time t'. A technical detail is that it can be apphed at time t' — e (which defines R^{z,t,t')) or at time t' + e 
(which defines R{z, t, t') and corresponds to the Ito prescription for the response, i.e 6 < x{t') > /6h{t') = 0). If quasi 
equilibrium hold (meaning the current J at ty^ vanishes, see Appendix B, the time translational invariant averaged 
functions Q{z, r) for large t' = tw and fixed r and R{z, t) should verify the exact differential relation: 



-drQ{z,T) = Td,R{z,T) (27) 

These relations can be generalized to discrete models as done in A. The response R'^{z,t) satisfies a slightly 
different equation: 



drQ{z,T) = -TdiQ{z,T)+Td,R+{z,T) (28) 



Wo determined both the distribution i?+(z, t) and R{z, r) numerically (as explained in B ). The response R'^{z, r) 
is plotted in Fig. 10. Note that for r = oo the above implies simply: 



R+{z,T) = d,Q{z,T) (29) 
and we checked that R'^{z, r) has the form <^ l/z^^f^^ expected from the above relation and (21). 
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R(z) 




FIG. 10. Plot of R+{z,t) [y coordinate) versus z for various t and comparison with the expected large t hmit l/z'"/^' 
("fit"). L = 200 and 5 10^ configurations. 



The Ito response R{z,t) satisfies (27) and is thus analogous to a probability current. Examining R(z,t) obtained 
from our simulation as a function of z confirms the above arguments based on the two well model: there is a (statistical) 
current flowing from the center region 2 ~ to the z > (and to the z < 0). In each disorder configuration R(z ~ 0, r) 
does not have to vanish (since a each local environment is not symmetric) but its average R{z = 0, r) must vanish by 
symmetry. We observe that in the simulation it is of the order 10~^ smaller than R{z = 1, r) which is consistent. 

Wc have plotted R{z = 5, r) in Fig. 11. It shows clearly a l/r decay. To be more accurate, we have plotted 
tIiitR{z = 5,t) in Fig. 12. It indicates that: 



R{z 



(30) 



with 6 « 1 (for a fixed z). This could be consistent with a decay of Q{z,t) ~ 1/ln^r for fixed z and a scaling 
^ In^ T in (27). However, Q and R are expected to consist of two parts with r dependent relative weights, one part 
scaling as z/ In r and a fixed part in z. Thus more work is needed to determine these functions more precisely. 

R(z=5,i; ) 




le+04 le+06 le+08 le+10 le+12 le+14 
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FIG. 11. Plot of R{z = 5, t) (y coordinate) as a function of r. Two simulations are indicated: L = 100 and 5 10^ configura- 
tions and L = 200 and 10^ configurations 
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FIG. 12. Plot oi T\og(T)R{z = 5,r) [y coordinate) as a function of r. L = 200 and 10^ configurations. 

Very recently (while this manuscript was in completion) a promising approach was developped by Cugliandolo, 
Dean and Kurchan (CDK)^" to obtain bounds which may permit to demonstrate that these quasi-equilibrium regimes 
exist in models such as Sinai's. Their approach is explained in Appendix C and some extensions and applications are 
given. For Sinai's model the bounds should be performed on disorder averaged quantities (since single environment 
ones do not converge). Doing this one can obtain two bounds of interest for Sinai's model: 

|< x^{t)>-< x{t)x{t')>-T [ R{t,t')\ < {<x^t) >)i/2 / ds(|^^|i/2 (31) 

Jt' Jt' ds 

and: 

\dt'Bit,t') + 2TR{t,t')\ < B{t,t'y/^\^^\'/^ (32) 

The first one was given in"**^ (though the problem of disorder averaging is not discussed) while the second is new. 
B{t,t') is defined in (10) and R{t,t') = 5< x{t) >/6h{t') = J dzzR{z,t,t') (with Ito's definition of the response). 

The function H{t') = J dx'P{x't' \00) (Tin P{x' ,t'\00) - U{x')) is the averaged free energy which satisfy an H 
theorem (it is always decreasing). The nice observation of"**^ is that since in Sinai's model H(t') ^ — hut' the above 
bound (31) imphes that the r.h.s bound for the integrated FDT violation is (Ini)^(i^/^ — t'^^^). Clearly ioi t = t — t' 
fixed it implies that (integrated) FDT holds, but it does even seem to imply that it holds further even for r ~ f'^/^. 

Note that the second bound (32) immediately implies that if there is a limit distribution lim^^.^oo Q{z,tw + T,tw) 
then it will depend only on t = t — t' (TTI) and the bound imply that FDT is verified (since dH{t')/dt' goes to zero). 
(Note that this is also the case for the Brownian motion). 

While these results are as inescapable as rigorous bounds, an explanation for this could be looked for in our previous 
arguments about periodic media. Indeed if our assumption (19) is correct, then one can choose In^ r <C <SC In^ 
and still have equilibrium and FDT (since the probability of having another accessible absolute minimum in the box 
Lq at distances 1 is vanishingly small (see )). Thus it is likely that in fact FDT will hold beyond what is shown 
by the bound, probably until r ~ with c < 1 (beyond that one enters the diffusion and aging regime, see next 
section). 

Note also for completeness the two single time bounds: 

\d,-^W)^\ < (33) 
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\d,< x{t'Y >\ < \< x{t'r >|V2|^^|l/2 (34) 

Thus we have found that an (equihbrium) TTI diffusion regime for the process z{t) persists within the quasi 
equiUbrium regime. 

One can also apply these bounds to the case with a an applied force / ^ /i, anticipating a little on the next Section 
III where the directed model will be discussed. Let us assume that H{t') ~ — < x{t) >~ — (i')^- Then (31) leads to: 

\<x^{t) > - < x{t)x{t') > -T f R{t, t')\ < t''{&+^^/^ - (t')(''+i)/2) (35) 

Expanding for r <C t' , one finds that the l.h.s must be smaller than t lt'^^~'^^^l'^ . Thus it shows the existence of an 
FDT regime for > in Sinai's model. The r.h.s goes strictly to zero as long as /i < 1/3 but if one actually divide 
the r.h.s by the expected scale of the l.h.s (^ t^'') one gets r/t'^^"'"'')/^. Thus a FDT regime (with X = V) seems to 
be possible for ^ < 1. 

This may seem surprising at first, since for /x > the system is driven. However it makes sense physically. Remember 
that for /i < 1 the particles will spend most (all for t^, oo) of their time in a well of release time t^, (see section 
III). There they have time to equilibrate. Thus in some sense the fact that there is aging (and broad distributions of 
release time) is intimately related to the fact that there can be equilibrium within a well. 

The above bounds also put constraints on the possible aging forms for a large class of models. This will be discussed 
in Section (IV). 



E. Aging and diffusion regime 



We will now present our results for P(n, t|n', i;')P(n', t'|no, to = 0) in the regime where t and t' are large and well 
separated. As stressed above the data is complicated to analyze. We have found numerical evidence for two regimes. 

(i) the diffusion regime 

The first one is the diffusion regime where z ~ (Inf)^ and Int' Int: 

^(^-(W^t(W'^] 

In that regime the three relative displacements x{t) — x{t'), x{t) — x(0) and x{t') — x{0) are of the same order of 
magnitude. 

(i) the aging regime 

Another regime was found by looking at the decay of Q{z 0, t, t') for z fixed and small (near ^ = 0) and t and t' 
large. By z fixed and small we mean in our numerical simulation n = n' but it means more generally that n is in a 

finite neighborhood of n' . For instance wc have foimd that the decay of Q{z = 0,t,t') could be as slow as desired by 
taking both t and t' to infinity. We have tried various dependences of the form: 



Q(z = 0,t,t^)^f[^] (37) 

The best fit was obtained for: 

h{t)r~.hit (38) 

This is illustrated in Fig. 13 where wc have used t,,, ^ t"' with < a < 1 a nd varied t with fixed a which is 
consistent with h{t) = Int, since a — In t^/ Int. For a = 1 one recovers Qeq{z = 0) as expected. Note that it is more 
difficult numerically to get correctly the small a behaviour (hence the curvature of the curves for small a on Fig. 13 
which is due to the smallncss of the available t^ ^ t"). We have also tried t,,, ^ t/2 in order to test a possible t/t^ 
dependence of this quantity (upper curve). This curve is in fact indistinguishable at large t from Qeq{z ~ 0). This 
strongly indicate that aging with a form t/t^ is unlikely. 

We want to emphasize that we could not rule out other forms more complicated than h{t) = Int. However, the 
form t/tw appears very unlikely for the measured quantity (e.g we haven't checked the behaviours of the moments). 
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FIG. 13. Plot of Q{z = 0,t,t°-) for various values of a which indicates that the aging scaling function is h{t) ~ Int. The 
curvature of the curves for small a at times not very large is due to the smallness of t". L = 200 and 5 10^ configurations. 



Wc have also computed Q{z, t, t"-) for several z fixed and large t. We find that the aging behaviour is also valid for 
z ^Q. In fact we find strong evidence for the behaviour: 

Q{z,t,t') = Qo{z)f\^^] (39) 

as is shown in Fig. 14 for < 2 < 10, where we have been able to collapse all the various curves Q{z,t,t') on the 
same curve Qo{z). Thus this curve should also be equal to Qo{z) = Qeq{z, r = oo) /Qeq{z = 0, r = oo). 

This behaviour seems to be consistent with a picture of aging in Sinai's model (at least for small z) where an 
equilibrated packet jumps as a whole out of its large well when t' ^ t"". This will be discussed further in Section IV. 
Though the collapse in Fig. 14 is perfect at small z (better than 10~^ in relative precision) it becomes poorer at larger 
z. This is probably due to the fact that for the available times one overlaps somehow with the diffusion regime. It is 
unlikely, but not ruled out, that this could be the sign of yet another regime (more work would be necessary). 



Q(z,t,tw=t °)/f(a)=Qeq(z) 
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FIG. 14. Plot of Q{z,t,t°-)/Q{z = 0,^,*") as a function of z for several values of a and for two times t = 10^° and t = 10^^. 
L = 250 and 10^ configurations. 

The fact that the temporal scaling is found to be the same h{t) ~ L{t) ~ Int in both this aging regime and the 
diffusion regime indicate that these two regimes may in fact be the same in this model. If that is indeed correct it 
does imply that the above two time diffusion scaling function (36) must have a singularity near the origin, e.g a delta 
function singularity. This remains to be investigated further. 

Note that in mean field, in short range models, one can also see that the diffusion scale and the aging scales actually 
coincide^^. Here it seems that the same happens, though for a different quantity Q{z ~ 0,t,t'). 

III. DIRECTED MODEL 

In this Section we study aging and diffusion in a simple model of directed diffusion with disorder where two time 
correlation functions as well as response functions can be computed simply. This model will appear as a particular 
case of a more general 'solvable model' which we will introduce in Section V 

In this model the particle can only jump to the right of the occupied site and quenched disorder is introduced by 
choosing at each site n an average waiting time r„ according to a given distribution P(t). The distribution of waiting 
times is further chosen with an algebraic tail at large times P(t) ^ C/t^~^^. This directed model was introduced and 
studied (for single time quantities) in^^. It was also shown'"''''^'^ that it is a large scale effective description of Sinai 
model in presence of a bias. The idea is that in presence of an applied force / the potential landscape {x, U{x)) is 
a biased random walk. There are thus some places where the walker goes back against the bias. This leads to rare 
barriers against the drift of size Eb with probability cxp(— /£"(,). Since the waiting time in these traps behaves as 
exp{—Eh/T) this is enough to generate dynamically an algebraic distribution of waiting times^^. Whether this model 
is also a good description of the biased Sinai model for two times quantities remain to be investigated in details. 

The model is defined by the Fokker-Planck (FP) equation: 

= H^mPmit) = Wn-lPn-l{t) " W„P„(0 (40) 

The rates at each site correspond to a mean waiting time t„ = 1/Wn- 

The Green function P{n,no\t — to) is defined as the solution of (40) with initial condition P{n,no\0) = Snno- Its 
Laplace Transform (LT) P{n,no\s) = Jq°° P{n,no\t)e~^* can be computed easily from (40). It reads: 

p(»,n.w = j^ n irk <"> 

fc— no 

Averaged single time averaged quantities are easily computed from: 



1 



P(0,0|s) = -— ^ = $(s) (42) 



P(n,0|.s) = $(.s)(l - .s$(.s))"-i (43) 
The long time limit is then obtained from the small s behaviour of the function $(s): 

$(s) ~ C7r/(sin/x7r)s''-^ < 1 (44) 

$(s) ~ 1/V + C7r/(sin^7r)s^-i 1 < < 2 (45) 

<!>{s)^l/V-Ds + CTr/{smnTr)s''-^ /x > 2 (46) 

This yields the several phases of the model: the subdiffusivc phase < < 1 where x ^ t^ (and zero velocity), the 
anomalous dispersion phase 1 < < 2 where there is a velocity V > Q but dispersion is anomalous D = co, and the 
diffusive phase /it > 2. We will be mostly interested here in the case < /U < 1. There the above result immediately 
yields a Levy diffusion front^^. 

We are now interested in the averaged probability that the particle advances by m between t' and t = t' + T: 
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Q{m, T, t') = {5{x{t) - x{t') -m)) = Y^ P{n + m, n, T)P{n, 0, t') (47) 

n>0 

This quantity is computed in the Appendix F. The result in Laplace variables is the following: 

$(.Sl) - $(.S2) 



Q{m = 0,81,82) 



{S2 - .Sl).S2$(s2) 



(82 - 81)82^(82) 

where si is associated to r and 82 to t'. In fact one notices that 

Q{m > 1, 8„82) = (— - p{m = 0, si, 52)) $(si)(l - si$(si))™-i (49) 
and thus the disorder averages factorize: 



Q(m > 1, T, f') = (1 - Q(0, r, t')) * P{m - 1, 0|t) (50) 

where the * denotes convolution over the variable t. This could be expected from the markovian and directed 
nature of the walk. The clock is set back to when they exit the trap. Thus for two times quantities all one needs 
to know is the two time averaged probability of not moving between tw and tw + t plus the averaged probability of 
diffusing by m during r. 

Let us first examine the probability of no motion. Explicit laplace inversion is simple on the asymptotic form for 
fjL <1 and it gives: 

Q{m = Q,Tt') = ^-^^^^^ f dxil-xY^'x^-^ ^F{T/t') /x < 1 (51) 

Jo 



where: 



One has: 



, sin(7rii) „ 

F{z) ~ — ^—^z-^" z^oo (53) 

/LtTT 

and 1 — F{z) ^ z^~^ at small z. 

It has a aging form as a function oir/t^. This is the manifestation of the Feigelman Vinokur trap model mechanism 
of aging. This expression is similar to the one obtained in^^ in an infinite range model. In the present case however 
we are also interested in the diffusion regime which we now analyze. 

We see on the above expression that if r is finite and ^ 00 the probability of being trapped is 1. There is no 
motion on finite time scales. Wc note that this directed model is not rich enough to contain information about the 
dynamics inside traps. The quasiequilibrium regime is thus degenerate Q{z,t) = SzfiSij). 

The above result also shows that at time t^+r the fraction of particules released by the well is of order (r/t^)^"'' 
and that the particles which are released do fast motion (as if the clock is then set back to when they exit the trap). 
In this model the later motion is not slower. Thus they will move typically by Sx ^ t^. 

Thus can easily estimate the moments: 



< {x{t^ + r) - x(t„))" > - (r/i^)i-''r"'' (54) 
Note that this gives, for the first moment: 

< {xit^ + r) - x{t^)) > - r/C ~ {ty, + - C (55) 
in the regime t << consistent with the known result for the first moment. 

We note that th ere is, in some sense an "aging regime" for each moment. Indeed one can impose that 
< {x{t) — a;(t'))" > is a fixed number while taking both — > 00 and r — > 00 provided: 
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This depends on the moment itself. This illustrates the strong non self averaging properties in one dimension. 

In addition to a self averaging aging regime which is confined to ^; = there is in this model a diffusion regime for 



z ^ and t ^ t' . It can be written as: 



Using the above exact relation: 



0(.,i,O = ^i=^[J,y] (57) 



Q{z, t, t') = Q{z ^ 0, t, t')P{z, t - t', 0) (58) 

where P{x,t,Q) is the single time diffusion front which takes a scaling form P{x,t,Q) ~ t~^P{xt~^) where P was 
determined in^^'^^ 
One finds: 

F[y, A] = (1 - X)->^F[^^]P[y{l - A)-^] (59) 

The response function can also be computed. Using the formula (AlO) of the Appendix A one sees that the matrix 
B = H and thus: 

and is thus independent of t'. In fact any response function can be obtained that way. For any observable depending 
only on time t: 

S < Ojt) >h _l d< 0{t) > 

Shf 2 dt ^ ' 

which comes from the directed nature of the model. 



IV. GENERAL DISCUSSION AND BARRIER MODELS 



In this Section wc first discuss aging and diffusion properties in low dimensional phase space in terms of scaling of 
barriers. To illustrate how these ideas work in practice, we study some cases which can be solved analytically, such as 
the directed model of the last Section with more general wide distributions. These ideas can in principle be applied 
to study aging in a wider class of models. 



A. aging, diffusion and scaling of barriers 

One way to formulate the question of aging is to ask what is the typical size of the next large barrier typically 
encountered by the particle after time tw (i.e at times t > tw)- The idea is that small barriers have already thermalized. 
One can usually consider that all barriers of size Tlntw — TC have already thermalized, where C is a constant which 
can bo chosen large but can be kept fixed as oo. (this will be true in 'reasonable' landscapes with fast enough 

growing barriers, e.g such that diffusion times are smaller than equilibrium times). Thus one must look at the next 
barrier encountered after which has a size bigger than Ef, > Tlntyj — C. For a relaxation behaviour of the type t/t^ 
to occur one must have that a typical particle in a typical environment (i.e that there is a finite fraction of particles 
and environments) must overcome a barrier within the range Tlntw — TC < Efj < Tlni^j + TC after time tw This 
is illustrated in Fig. 15. If typically there is no such barrier, e.g the next encountered barrier is always larger, then 
the particle is typically thermalized iov t/tw = c and thus one expects: 

lim Q{z, ctw, tw) = Qeq{z,T = oo) (62) 

which seems to hold for the Sinai model from our simulation results (see section II). This situation is illustrated in 
Fig. 16. 
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On the other hand, if the particles typically encounter such barriers, then one can have a scaling as t/tw For that 
one needs an appropriate degeneracy of the barriers. This is the case for the directed model which has en exponential 
distribution of barriers e""^. At there is a finite probability <^ Q-a{Eb-Tintw) ^j^g^^ next barrier is Eh- A more 
precise analysis is presented in the next section. 



ln(t„) 




FIG. 15. Barrier condition which leads to aging as t/tw 

To have aging in t/t^ there is of course no need to have a scenario a la Feigelman-Vinokur and also no need for 
randomness. In the next section we construct a simple model which exhibits a variety of aging and diffusion regimes, 
but which satisfies the above mentionned property. 



ln(t„) 




FIG. 16. Barrier in Sinai's model 

The property (62) does not mean that there is no aging, but rather that it is degenerate in t/tw In fact (62) suggest 
that there must be a scaling which gives some aging. One can always define a function t = gdtw) such that: 



Q{z = 0, gdtw), tw) =c < c < Qeq{z = 0, r = oo) 



(63) 



t'^'^^ and 



The question is how does the function qcit,,.) behaves for large tw A natural possibility is that gdtw) 
indeed we now argue that t his is the cast' in Sinai' s model. Note that more generally one can define for each z a 
fmiction g{c, z, tw) such that Q{z = 0, g{c, z, tw), tw) = c. If the large tw behaviour of this function is the same for all 
z then one has an aging regime with a non trivial z dependence as in the models of Section (V). 

In the symmetric Sinai model the question of the next large barrier to be encountered is difficult to answer analyt- 
ically but the simple following arguments can be made. Typically at time t^ the particle is within a valley (see figure 
16). Let us call Umin the lowest energy point of the valley. The walls of the valley are at least Umin +T\ntyj — TC. 
The particle sits at a point which is within U < Umin + TC in energy (since the valley is typically thermalized, sitting 
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any higher has a vanishingly small probability). 

Geometrical construction of possible such valleys in the energy random- walk landscape of Sinai's model shows that 
typically the next barrier, i.e the size of the smallest wall is of the order of (1 + a)Tliitw where a is a fluctuating 
positive random variable. This can be seen from the fractal nature of the landscape. This suggests that a dependence 
of type \nt/ In tyj, i.e: 

lim Q{0,t,t^) = F[^] (64) 

which is indeed observed. Note that one also expects, from the argument after (63) that an aging scaling form 
should exist for any finite z: 

lim Q(z,t,t^)=F[z,^] (65) 

Thus to predict aging properties one needs to know the statistical properties of sequence of barriers effectively 
encountered by a particle. It would be appealing to relate it to the geometry of the environment. A first step would 
be to define in a one dimensional landscape the sequence of next largest barrier encountered starting from an initial 
point. Let us call them G„ = Eb{n). Two examples of landscapes are shown in Fig. 17 and Fig. 18. Fig. 17 represents 
e.g the directed model of last section, with an arbitrary distribution P{E) of barriers (of height Ei corresponding to 
waiting times = eKp{—Ei/T)) which generalizes the algebraic waiting times distribution. The important question 
is how does the series G„ = Ef,{n) grow typically with n. This will already give a rough idea of the type of aging 
one can expect. In particular it may determine the aging scaling function h{t). In Sinai model the series follows a 
(random) geometric progression Eh{n) ^ exp(cn) as illustrated in Fig. 18 which leads to the observed behaviour (64). 
In order to have at/t^ behaviour one needs that Ei,{n) grows typically as cn which is what happens in the directed 
model. Note that the progression of barriers Ei,{n) with n depends only on the topology of the lansdcape and for 
instance of distances along x (different landscapes can have the same Eb{n)). 




FIG. 17. Sequence of next largest barrier in simple models (see text) 




FIG. 18. Sequence of next largest barrier in Sinai's model starting form point 

Of course to predict more detailed properties one may need to know more: i.e the sequence of barriers effectively 
encountered by the particle. For this one needs also to know about the well depths. Let us consider e.g. a non directed 
model but with a landscape as in Fig. 17 where all the wells have the same depth. There the packet will thermalize 
in a large region (bounded by the next largest barriers) of size L{t). This will result in a geometrical prefactor e.g 
Q{z = 0,t,t') ~ L{t)~^ F[h{t) / h{t')]. If the well depths grow sufficiently fast then the packet will be concentrated at 
the bottom and then one may expect a reduced, or even finite prefactor L{t) = est (as in Sinai's model). 
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Thus if one knows the statistical properties of the sequence of next largest harrier one knows a lot about the aging 
form. In fact if these barriers grow sufficiently fast with n (faster than n) we expect that this information is enough to 
determine these aging form entirely (provided the wells grow also fast enough - see remarks of previous section). This 
suggests, e.g in Sinai's model, a program to study aging exactly (though it is technically difficult). Let us illustrate 
these consideration on models where this construction can be done easily. 

classification of aging forms and constraints from the CDK hounds 

Before we do so let us present some general considerations about the aging functions h{t). Remember than there 
is some gauge freedom in choosing them since aging quantities arc determined by a fixed ratio h{t)/h{tw) = c. Thus 
e.g the choice h{t) ~ Hs as good as h{t) = (with the proper change in c). 

An obvious classification is to distinguish between three classes of functions h{t). 

(i) h{t) grows slower than t (or any power of t). This corresponds to fast growing barriers. The condition 
h{t) / h(tw) = c can also be expressed asr — t~tw^t^ H{c, tw). Examples are: 

/i(0~lnt tr~.Tr~.tl, (c> 1) (66) 

/i(i) ~ e(^"*)° f ~r~i^exp(^(lnt^)^-'*) 0<o<l (67) 

/i(i) ~ e('" i ~r ~i^(lni^)'"'^ (68) 

(ii) simplest aging behaviour h(t) = t. 

(iii) h{t) grows faster than t (subaging). Then one must have t = t — tw ^ t^, such as in the following examples: 

/i(t) ^ eC"*)" r=— — ^ a>l (69) 

^ a (lni„)''-i ^ ' 

h{t) ~ e*° T ~ — ti,-'' < a < 1 (70) 

(71) 

The CDK bounds discussed in Section II also put general constraints on possible aging forms assuming they have, 
as in mean field, a non trivial < X ^ 1 < oo (see Appendix B for definitions). Assume indeed that as in mean field 
all correlation functions and X are functions of h{t')/h{t). Then bound (32) implies that: 

dln/i(t') ^^|^^|i/2 (72) 



dt' ' dt 



Thus if H{t') ~ i' " one has that 



h{t) < exp(t(i-")/2) (73) 
which excludes a large class of subaging behaviours (Sinai corresponds to a = 0). 



B. sequence of largest barriers 

It is simple in some cases to determine exactly the distribution of the sequence of next largest barriers. 

Let us look again at the directed model defined by a set of successive barriers Ei, identically distributed with a 
distribution P{E). We denote H{E) = /+°° P{E')dE' = Proba(£;' > E) 

Let us estimate the probability density that the sequence of successive next largest barriers (see Fig. 17) be 
Eq,Gi,G2, --Gn- It is by definition: 

Q(i;o,Gi, ..G„) = ■■ Proba(i;i < So, -Eu, < Eo,Ek,+i = G2 > Eo,Ek,+2 < Gi, ..Ek,+i = G2 > Gi, ... (74) 

fei— 0,00 fc^— 0,00 

■■Ek„+i = Gn> G„_i) (75) 

This yields immediately: 

^ . , p{G,)e{G, - Eq) piG2)9{G2 - G^) p(G„)g(G„ - G„-i) 
Q{Eo,..Gr.)=p{Eo) jj^^ jj^^ ... jj^^^ (76) 
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It is easy to see that the exponential distribution p{E) = exp{—E)0{E) has a special property. Indeed, in that case: 
Q{Eo, ..Gn) = e{Eo)e{Gi - Eo)e{G2 - Gi)..6{Gn - G„_i)e-^" (77) 
but this can be written as the product: 

Q{Eo,..Gn)= n MGi-Gi_i) (78) 

1=0,1, ..n 

with Go = Eq and G_i = 0. This shows that the difference between successive largest barriers have the same 
exponential distribution p{E) and are independently distributed. Thus: 

i=0,n 

The central limit theorem can then be used and leads to a Gaussian distribution for the variable (G„ — n)/^/n. 
This remarkable property of the exponential distribution allows to understand the Feigelman Vinokur scenario for 
aging as a linear (random) growth of the next largest barrier Ef,{n) = G„ ~ an 

This remarkable property of the exponential distribution also allows to obtain the general solution of the problem 
of finding the probability of a sequence of largest barrier. Defining ^{E) = - In dE'P{E') one has: 

P{E)dE = de-^^^^e{E) = ^'{E)e-^^^^e{E)dE (80) 

and thus (76) can be put in the form: 

Q(Go,..G„)dGo..dG„= (81) 
d$(Go)...d$(G„)^?(Go)^(Gi - Go)0(G2 - Gi)..0(G„ - G'„_i)e-«°"-°''-i)+(°"-i-^"-=)+-+(^i-^°)+^''' (82) 

and the sequence of variables $(Go), ..$(G„) can be constructed as a "random walk": 

n 
i=0 

where the wi arc a set of independent variables, identically distributed with an exponential distribution P{w)dw = 
e~'^9{w)dw. There are several consequences. First the distribution of each $(G„) is thus a Poisson process: 

Q(G„)rf$(G„) = ii^e-*(«") (84) 
n! 

and one has the central limit theorem for large n: 

^(G„) - n 

where w is a centered gaussian random variable of unit variance. 
One can apply these results to various cases: 
(i) algebraically growing barriers: 

P{E)dE = aE"-'^e-^''e{E)dE which corresponds to <P{E) = E". Then one has: 

,1/a 



(85) 



G„ ~ n^/" (86) 

Note that the randomness (disorder) is only a subleading correction, 
(ii) exponentially growing barriers: 

P{E)dE = (1 + E)-^"-+^^e{E)dE which corresponds to ^{E) = alnE. Then one has: 

G„ ~ e"/" (87) 

The variable G„ has a log normal distribution. 
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C. consequences: aging properties 



Having determined exactly the sequence of next largest barriers one can construct a toy model which will mimic 
the exact diffusion process. Wc claim that if barriers grow fast enough (faster than £"„ ^ n) it will give the exact 
aging behaviour. It can be applied to a variety of landscape, but let us apply it here only to the previously considered 
directed model. 

The toy model amounts to approximate the quantity Q{z, r, tw) at 2; = (i.e the probability to remain still between 
tu) and tw + t) as: 

+00 

Q(t, t^) = Q(0,i^ + r, t^) = e{Gk-i <T\nt^< Gk) exp(-re-^'=/'^) (88) 

fe=0 

This means that at all barriers smaller than T In have been overcome and that the only relaxation process is 
to go over the next highest barrier. This is illustrated in Fig. 17 which is adequate for the directed model. It thus 
supposes that the probability to be in another well is zero, which is correct for the directed model if barriers grow 
fast enough (which means distribution of waiting times wider than power laws) 

This can be rewritten as: 

+00 

Q{T,t^) = Y,i<^{^r. - *fe-i) - e{^^ - $fe))exp(-Te-+*"'(*'=)) (89) 

fe=0 

where $fc = $(Gfc) and = $(rint^). Because of the above statistical properties of the sequence this yields 
upon averaging: 

dcbY^QkW rfw;e-"'(^($^-<^ + «;)-e($^-(/)))exp(-Te-+*"'(^)) (90) 

fe=o 

and we can use (84) namely that J^t^ Qk{4') = 1- In the large time limit one can shift the integrand globally by 
+ $w without edge effects (though one must be careful not to shift terms independently because of divergent 
integrals). Thus: 

Q{T,t^)= d4> du;e-'"(6'(-(/. + u;) -6'(-(/>))exp(-Te- + *" ('^+*"')) (91) 



This can be rewritten, after integration by parts as: 



Q{T,t^)= / di;e-"exp(-Te-** (-+*-)) (92) 
Jo 

This is our general result for this toy model. 

Let us estimate (92) first for models with fast enough growing barriers, i.e faster than $(i?) ~ E. We introduce 
the aging scaling function h{t) = g^C^'"*) and thus $(a;) = lnh{e^^-^) and i>~^(y) = Tln/i^^(e^). One can rewrite: 



Q{T,ty,)= rf«;e-"'exp(-e-'?'^('"'^'*»)) (93) 
Jo 

with S{w, T, tyj) — ^~^{w+^{Tlntyj)) — T\nT. It turns out that if barriers grow fast enough the function exp(— e~'^) 
acts exactly as a theta function 6{S) and the result is simply: 

Q{T = t-t^,t^) = (r » t^) (94) 

This can be shown by a careful examination of the asymptotics. The point is that if one takes r and tw to +00 such 
that the ratio h{tyj)/h{T) is fixed then the result is (94). Since barriers grow fast enough one has \imh(tw) / h(t) = 
Y\mh{tu))/h{T). A more correct version of the above statement is that: 



lim Q{t = t -tw.tw) = y (95) 
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It is interesting to note that one has exactly: 



h{tw) 



(96) 



and thus the above resuh follows from the exponential distribution of w. Note that this is like taking the T = 
limit and thus if barriers grow fast enough we are dealing with a T = fixed point. Note that our result can also be 
rewritten quite generally as: 

Q{T = t-tM = {Voo^i^if (97) 
it„ P{r)dT 

as a function of the waiting time distribution r = e^/-^. In particular we have determined that the aging function 
is exactly: 

^^'^ = f+°° if w ^^^^ 

Jj P{T)dT 

One can also define the probability that at time tw the particle is next to a barrier (i.e in a site) of waiting time f, 
i.e: 

r+oo 

Q{T,t^)= / dfPt„(f)e-"/^ (99) 
Jo 

(100) 

Then the above formulae yields that this "aged" waiting time distribution is: 

PtAf)df = ^t^^(Mr) - h{t^)) (101) 



This also takes the simple form: 



Jt,.. P{T-')dT' 



which we have shown holds exactly. 

We can now apply these results it to several cases. 

(i) algebraically growing barriers: 

$(£) = E"- with a < 1. Then one has h{t) = exp((Tlni)«). 



Q{T = t- t^,t^) = exp(-[(rinf)'^ - (Tlni^)'^]) (103) 
(ii) exponentially growing barriers: 

P{E)dE = (1 + E)-^"+^^e{E)dE which corresponds to ^E) = a\nE and h{t) = {Tint)". Then one has: 

Q{r = t-t^,t^) = i^^r (104) 

Let us now apply the above model in the case of linearly growing barriers {^{E) = E). There we know that the 
assumption that the whole packet is concentrated nest to the next highest barrier is certainly not exact. Indeed since 
there are many barriers within Tlnt^ — C < E < Tlntyj + C the packet will be spread out (barrier degeneracy). It 
is instructive though, to see how well this toy model does in that case. 

Setting y = r/tw one gets: 

/■+00 pW 

Q{y)= / dwe-"^ / #exp(-ye-'^/^) (105) 
Jo Jo 

After integration by part and change of variables: 
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Q{y)= dw;e-«'exp(-ye-"'/^)=Ty-^ / ^A^e"^ (106) 
Jo Jo A 

For large y it behaves as: 



Q{y) ~ TT[T]y-^ (107) 

which is to be compared with (53) (T is exactly fx, T = /j,. from the relation between waiting times and barrier 
heights). The toy model gives exactly the exponent, but not the prefactor. We note however that at small T the 
prefactor becomes exact which is in agreement with the fact that for faster growing barriers the toy model becomes 
exact. At small y the toy model however docs not yield the non analyticity at t << tw, which thus entirely originates 
from the initial (fractal) dispersion of the packet at tyj over several wells. 

This confirms our physical picture: only maxima start playing a role when barriers G„ grow faster than n. The 
case Gn ~ n is the marginal case when finer information is important. 

Note that the exact result (51) can be put under the form: 

r+oo 

Q{T,t^)= rfu(;6(u)e-"^/*» (108) 
Jo 

, sin(7rT)e-« r dt ^ t 



to be compared with the toy model which has a <p{u) = Tu^~^6{0 < u < 1). Thus one can see that for large w 
the is correct, however the barriers which are not large, or the one at w < (smaller than tyj) start playing an 
important role. 

The role of smaller barriers can be illustrated as follows. The correct distribution of the waiting time f of the site 
where the particle is at is computed in the Appendix. It reads: 

(f)rff = /'e-^(^-«)^ (110) 

TT tyj T Jq T\}l\ 

This should be compared with what the barrier model which gives: 

PtAr)dT = „-C-^f+^^9{f-ty,) (111) 

Zy) T 

from the above fomnlae (102) and the identification f = e^l"^ . As above, the toy model gives the correct power 
law dependence for large waiting times (the prefactor itself becoming exact when ^ 0). There is however an 
accumulation of smaller barriers effectively seen by the particle, which is not captured. 

Let us conclude this section by noting that this type of distribution of next largest barriers can be used to analyze 
a large variety of models. For instance the analysis of the symmetric waiting time model in d = 1 or of directed model 
with several branches will be quite similar but goes beyond this paper. 



V. A SOLVABLE MODEL WITH AGING AND DIFFUSION 



We will now present a solvable model which exhibits simultaneously a non trivial aging regime and a non trivial 
diffusion regime. The third regime (the quasi-equilibrium one) is degenerate (it is reduced to a point). More properly 
this is really a class of solvable models, and we will only study a few. 

In view of the discussion about barriers of the preceding section the best way to construct a one dimensional 

diffusion model with a non trivial aging regime (i.e such that there is a finite probability that z ~ x{t) — x{t^) remains 
finite when both t and tw are large) is to make sure that the next largest barrier seen after time tw is equal to the 
previous one plus a constant (for aging as t/ tyj). A natural landscape is thus to look at a succession of barriers 
E" ^ n. However one wants the valleys also to become deeper and deeper, otherwise the thermal packet will be too 
extended. Thus a natural choice is also to suppose valleys to scale as ~ — n, a landscape represented in Fig. 19 

It turns out that the continuous version of this model, as well as some generalizations, can be solved exactly in a very 
simple way. 
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FIG. 19. Aging model studied in the text in the continuum Umit 



A. the general model and its solution 

Let us consider the following one dimensional diffusion equation: 

dtPix, t) - d^{D{x){d^P{x, t) - F{x)P{x, t))) (112) 

This model corresponds to diffusion in a landscape with both barriers Ei){x) and a potential U{x) (valleys) such 
that: 

D{x) = e-^--^^) U'{x) = -F{x) (113) 

here and in what follows we will set the temperature T = 1 for convenience but it can be easily put back in. The 
potential is defined by the fact that the equilibrium zero current measure is Peq{x) = exp(— [/ (x)). The general model 
(112) cannot be solved but there is a particular case which can be solved easily, which is: 

U{x) = -Et,{x)/2 (114) 

The model is defined by giving a function ^(x) such that: 

^U{x) ^ ^-EUx)/2 ^ ^^^^ ^^^5^ 

Then the diffusion equation becomes: 

dtPix,t) = d^mx)d,{^{x)P{x,t))) (116) 
and is easily solved since one can define a new variable u and a new probability G{u, t) such that: 

^ = P(x,t)dx = G(u(x),t)du(x) (117) 

ax 9{x) 

In the variable u the system is a free diffusion problem ! 

dtG{u,t) = dlG{u,t) (118) 

Clearly for the choice $(x) = and thus Ei,[x) = 2x this model is some continuum limit of the one represented 
in Fig. 19. It will have the same expected large time properties. Note that the quasi-equilibrium regime will be lost 
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in this limit (i.c reduced to a delta function, see below) as the size of each well will become infinitesimal. It would be 
nice to be able to solve directly the model of Fig. 19. to obtain also the small times behaviour. In any case the results 
presented here for large times will be the same (they arc not an artefact of the continuous limit). 

One must define carefully the boundary conditions. There are basically two choices that we will study: 

free boundary conditions 

One can study a barrier landscape Eb{x) defined from — oo < a; < oo. This could either be a random landscape 
U{x) or a deterministic one with growing barriers in d = 1 which we then should choose symmetric for definiteness 
for instance U{x) = —\x\. Then a natural definition of m is : 

u{x)= r dx'e-^^'''^ (119) 



and — oo < m < oo. 

The Green function is thus simply unbounded diffusion 



1 (u-un)^ 

G{u, t\uo, 0) = -^=e — (120) 



which yields the Green function for the original equation: 



P{x, t\xo, 0) = -^e-^^^'^e-— — (121) 

We will be interested in the probability of displacements z between t' and t: 

Q{z,t,t')= j dxP{x + z,t\x,t')P{x,t'\xo = 0,0) (122) 

which formally reads: 

Jo J_x 

'4^^^/t'{t-t') 
Another convenient form is to write: 



Q(z,t,t')= dx g-c/(x)-[/(x+^)g 2 4(,-,') (123) 



duidu2 S{z - (x(ui + U2) - x(ui)))e~^~^(^ (124) 

4TT^t'{t - t') 

where x{u) is the function implicitly defined by (119). 
reflecting boundary at u = 

Since we will be interested in landscapes with growing barriers such as depicted in Fig. 19. it is useful in certain 
cases to introduce a reflecting boundary on the left. One can choose for definiteness a landscape such that U{x 
—00) +00, define u{x) = dx'e~^^^ ^ and use a reflecting boundary at u = 0. This naturally avoids the 
particule being either in one half space or the other. For instance for the lanscapc E}y{x) = 2x, for which u{x) = 
the reflecting boundary is at x = —00. Because of reflecting boundaries one must choose the free propagator with 
dy,G{u = 0) = 0, i.e: 

G{u,t\u',t') = , +e"^^^^) (125) 

^Awit - 1') 

We will be interested in: 

Q{z,t,t') = j dxP{x + z,t\x,t')P{x,t'\xo = -oo,0) (126) 

where we have chosen for convenience the initial condition at = — 00 (and thus uq = 0). This is a purely 

technical point on the definition of the model and has no bearing on the physics. Indeed since we have chosen 
Eh^x —00) —00 and the initial condition will be immaterial since it takes only a finite time in this model to 
reach finite x values and we are interested only in the later (long time) behaviour. Thus using: 
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Using P{x,t\x',t') = ^{x)-'^G{u{x),t\u{x'),t') one finds: 

Q{z,t,t') = / dx g-c/(x)-[/(x+^)g-H^(e J(l^F^ +e ) (127) 

J 4:^T^/t'{t-t') 

Let us conclude this section by indicating that a more general case can be solved (see also D), i.e brought back to: 

dtG{u, t) = dlG{u, t) - vduG{u, t) (128) 

which corresponds to: 

U{x) = -{Ei,{x)/2 + v r dye'^'^y^/^) (129) 
Jo 

These models corresponds to either valleys and barriers scaling differently and we will not study these models here. 
We note that the directed model of section III can be seen as a particular case of the class of models introduced here. 
Indeed the case (128) corresponds to the diffusion in x given by: 

vt' + VVw = u{x) - u{xo) = I e-^^"^ Ux' (130) 

J xa 

where w is a normalized gaussian variable. In the fully directed case (large v) and for a judicious choice of the 
(random) U{x) one can recover the directed models. 



B. solution for linearly growing barriers 

It corresponds to the model of Fig. 15. Again one can either consider the symmetric landscape taking U{x) = — \x\ 
and thus x{ti) = sgn{u) ln(l + |u|). Or one can take the "half landscape" with a reflexive barrier at xq = — oo which 
we will discuss first. 

The above formulae give in that case, for the single time packet: 

P(x, t\xo,Q) = -^=e^e-^ (131) 



The generating function of moments is < >= T[{1 + A)/2](4t)-*''/^/v^. Thus one has: 

<x{t)>=^\n{t/c) (132) 
<x^{t)>- <x{t)>^^TT^/8 (133) 



X 



where c = ln7 (7 is Euler's constant). In fact the packet has a limit shape as one can see by performing the shift 
i lii(4f) + X. Then the distribution of x is asymptotically time independent: 
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P{x,t) r-. ^e^e-^^^ (134) 

Thus the packet has a constant width and simply spreads over a few wells (in Fig. 15) with its center moving 
logarithmically towards the right. Thus, since there is some degeneracy of barriers one expects a.t/tw aging behaviour. 
Indeed one gets for the two time packet: 

Q{z,t,t')= / duu e~SF(e *(*-*') + e *(*-*') ) (135) 

This yields to an aging form for the distribution of displacements between t' and t, given by: 

which, in the variable w = , is the sum of two Lorentzian of width F the first one centered around w = \ and 
the second one which is its mirror image (with the mirror at w = 0). This result is natural considering that the ratio 



29 



of two independent gaussian variables with unit variance v = v\/v2 is the Lorentzian P{v) = l/{Tr{v^ + 1)) and that 
one has: 

z = x{ui + U2) - x{ui) = ln(l + — ) = ln(l + Tv) (137) 

where wi = VFvi is the positive gaussian variable representing the diffusion process on the half line and U2 = \/tv2 
another gaussian variable representing the later diffusion process, constrained so that the sum ui +U2 remains positive, 
hence the mirror. The total weight in the second packet is p = (arctanr)/7r and in 1 — p the first. 

Thus the first Lorentzian packet in (136) corresponds to particles which have remained in the region x(t') while the 
others have crossed to the mirror at a; = —00 at least once and came back to that region. We are mostly interested 
with this first packet (which always contains a fraction > 1/2 of particles) but the second packet will always be 
there as a mainly technical feature of the model. If wc chose instead the symmetric environment, these would be two 
separate aging packets one around x{t') and the other around —x{t'). Aging then occurs only within each packet 
(there are diffusion events between packets) while here one has aging in all the packet. 

Indeed the above distribution (136) clearly exhibits an aging scaling form of the type: 

Q{z,t,t') = Q{z,^) = Qiz,^^) (138) 

thus h{t) =tm this model. The function Q indeed depends on the times only through T which itself can be written 
as: 

r = /[|Il (139) 

where the form of f{x) and h{t) is univoquely determined as h{t) = t and f{x) = Vl — x. Thus, as in mean field^^ 
there is a singularity at the beginning of the aging regime: the /? exponent is equal to /3 = 1/2. 

It is interesting to note though that at the beginning of the aging regime F << 1, there is anomalous behaviour of 
the moments z", because the Lorentzian has diverging moments: 



/+00 J 
(ln(l + rw))"- 
1/r 1 



(140) 



One finds in particular < z >= l/21n(l + T^) = l/2\n{t/t') exactly. 



<iz-<z >)" == / duu"( — + — ) (141) 

2TrVT+T^J-oo coshM-^jipj cosh[u] + ^^' 

One finds that for n > 1 < (z— < z >)" rr[n + 1] (i.e P{u) = Fe"") which is strong intermittence. 
For widely separated time scales there is also a diffusion regime in this model. At large t one gets: 

z = ln(i;) + ln(t) - In(t') (142) 

where w has a Lorentzian distribution. Thus the diffusing packet has a finite size (note that since for a Lorentzian 
< In(ti) >= (inversion symmetry) one recovers the above result). Thus there is a diffusion regime defined by 
z,t,t' 00 with z/\nt and Int/lnt' fixed. It has no thermal fluctuations (it is completely determinist) and the 
diffusion scaling function reads: 

g(.,,,<V.~^*(jij-(l-5^„ (143) 

As in SR mean field models the aging regime smoothly merges in the diffusion regime^^ 

This model thus contains both an aging regime in t/tw and a diffusion regime with a different scaling. Only the 
FDT regime cannot be seen, since it has disappeared in the continuum limit. 

Finally note that this model seems to violate the quasistatic assumption. Indeed in a box of finite size L, or with 
periodic boundary conditions and taking ^ 00, the packet will converge towards the equilibrium measure. Since 
U{x) = —X this will result in an equilibrium packet oi finite size. However the dynamical regime t — t' t' leads only 
to a delta funtion packet at 6{z), thus totally different from the equilibrium one. We expect this feature to persist for 
the model of Fig. 15, i.e it is not an artefact of the continuous limit ! 
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Other quantities can be computed in this model. Let us give some examples: 



calculation of separation of replicas 

One can also compute the separation of two thermal replicas which are allowed to split at time tyj. This quantity 



was studied in 
One has: 
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and one finds: 



Q2{z,t^,,T) = j dydxP{y,t\x,tw)P{jJ + z,t\x,tw)P{x,tw\xo = -oo,0) (144) 

r + co r + co 

= I dui / duuiG{ui,T\u)G{uie'',T\u)G{u,tw\0) (145) 
Jo Jo 

G(«,t|«') = (exp(-^^-^) + exp(-^^-^)) (146) 



Qi{z, tw, r) = —{-r-z ^ r2 + ( ^^2^r2 ) (1^7) 

TT (e^ — m)^ + {e^ -\- my 



V=ATpl±A (148) 

with T = t — t' . Note that this distribution is symmetric under z ^ —z because of the relation + = 1. For 
large t ^ t^ one finds that the distribution of a; = goes to a fixed half Cauchy unit distribution 0(a;)l/(7r(l + x"^)). 
Thus for large time we find that the two replicas evolve within a finite distance, but this distance is larger than the 
dynamical qea which is zero for this problem. 

Note that when r/tu, <C 1 one recovers exactly the previous Q{z, t, t'). 

non trivial FDT violation ratio ? 

It is interesting to know if one can find finite dimensional models with a non trivial FDT violation ratio X as in 
mean field. 

Here one can also compute the response to an additional field, i.e the remanent magnetization decay. The calculation 
is indicated in the Appendix E. Though it does appear that in some sense this model has a non trivial FDT violation 
ratio X similar to mean field we were not able to exhibit it in a clear way. If one looks at the finite fraction of the 
packet which has not touched the refiecting boundary, it has clearly a non trivial X. But on any global quantity we 
have looked at the boundary effects always introduce a cutoff which changes the expected result. We do not know 
if this is a purely technical limitation, and if the model can be improved to really exhibit a non trivial and properly 
defined X or if this is a more fundamental limitation. We still present some of the calculations in the Appendix E for 
the brave who is encourage to improve on it. 



C. solution of more general deterministic model 

1. barriers growing faster than linear 

One can study cases where Eb{x) ~ with & > 1. From the Section IV we expect that aging should have simple 
properties. For the "half landscape" model, using u{x) ~ exp(a;'') for x > and a left reflecting wall at m = 1, one 
finds an aging regime with no thermal fluctuation: 

z = {\u{\fevi + s/t^ViY'^ - (In(Vt'ui))^/* (149) 
z = (In ^^tfl^ - (In sft'fl^ (150) 

where v\ and U2 are uncorrelated normalized gaussian variables. We have performed an expansion, in the regime 
of interest t^t' e.g: 

(ln(VF^i))i/''^(lnVf)^/''+ ^'ti-''^^ , (151) 
\ \ ^)) \ I (inVF)i-Vf' ^ ' 
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and dropped the contributions of the noise parts vi and V2 which vanish in the Umit of large t and t' since b > 1. 
Thus aging becomes thermally deterministic as: 

with h{t) = exp(C(lni)^/'') and C = (1/2)^/''. These results are strikingly similar with the result (96) of the 
previous Section. In fact the result (152) is also similar to what was found in mean field in^^. 

There is also a diffusion regime, which is identical to the end of the aging one (it merges smoothly into it). Indeed: 

. = (lnv^)V''(l_(^)i/^ (153) 



2. barriers growing slower than linear (suhaging) 

Similarly, in the subaging case where Eh{x) ~ x*" with 6 < 1 one finds: 

z = ln(l + Vtv2 + Vfviy^'' - ln(l + Vfviy/'' (154) 
Performing expansions it gives: 

^ = 1 ^ (155) 

and thus 

z=\ (156) 

where v is again a variable with a Cauchy distribution (we have not determined the exact form using the wall). 
This is compatible with a Lorentzian aging packet but with: 

r = iV^ (157) 



consistent with: 



^^^''W)^ Kt)=eMC{\r,tf") (158) 



with a < 1, again performing an expansion in small r/tw and with /? = 1/2. 



VI. CONCLUSION 



In this paper we have investigated two time quantities in several one dimensional diffusion models with random 

and non random environments. These quantities are functions of the waiting time t' = t^j, after the initial localized 
condition aX t = and a later time t. The (averaged) distribution Q{z, t, t') of relative displacements z = x{t) — x{t') 
between t and t' was studied. Part of this study was numerical (in Sinai's model) and we have reached times up to 
10^^. Our results showed that the times reached in a previous simulation'''* were vastly insufficient. Our conclusions 
are different from the one of''^. Part of the study was analytical: we have computed two time quantities for a directed 
model related to the biased Sinai model and we have introduced several new models which can be studied analytically. 
Our main results are the following. 
We have identified three generic regimes for large times tw ^ 00, t ^ 00: 

(i) At small separations t = t — tw & quasi-equilibrium regime. Evidence for that regime was found in Sinai's 

model. In that regime Q{z, + r, t^) reaches a limit Q{z, t) for t^, 00. We have argued, and checked numerically, 
that this distribution has some peculiarities. For large r it does admit a limit Q{z) but this limit exhibits an algebraic 
tail originating from rare configurations of the disorder. The moments of the relative displacement < z"{t) > with 
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n > 1/2 grow unboundedly with r. We have also proposed an expression for the distribution Q{z,t) based on 
arguments on periodic media. 

We have also concluded, from our simulations and from physical arguments, that in Sinai's model usual equilibrium 
theorems hold in this quasi equilibrium regime (TTI and FDT). We have obtained a generalized expression of these 
theorems to probability distributions such as Q{z,t,t'), and shown that in this regime Q{z,t) and the response 
function R{z,t) obey and exact differential relation. That these theorems should hold in that regime is confirmed by 
recently obtained rigorous bounds, as we have discussed. This unveils an interesting situation of a quasi equilibrium 
regime with a lot of internal structure, wide fluctuations, and internal logarithmic diffusion (moments growing with 
r) which calls for further studies. 

(ii) at large time separations, L{t) ~ L{tyj) there is a diffusion regime. There the displacements scale as x{t) ~ 
x{tw) ~ L{t) and there are scaling forms for the probability distributions. In the model of Section (V) we have 
obtained this regime analytically. 

(ii) finally there is an intermediate aging regime. One should first look at this regime in the probability of staying 
in a finite neighborhood z of the same point between t and t' which is generically of the form 

Q(z,t,t') = F[z,^^] (159) 

In Sinai's model with a bias (and in the directed model with algebraic distribution of waiting times) one has aging 
with h{t) = t. If the waiting times are even more widely distributed, we find (159) (for z — Q) with a large class of 
functions h{t) ^ t. Similarly we also find this behaviour in a solvable model in (V) where a large class of functions 
h{t) (including subaging h{t) <C t) can be obtained. 

In the symmetric Sinai model (without a bias) we have found strong numerical evidence for the aging behaviour 
(159) with h(t) ~ Ini (for small finite z). In that model we have even found a more striking result: 

Q[z,t,t') = Qo{z)f[^] (160) 

i.e a decoupled form for the aging regime. This suggests an interpretation of the aging regime in Sinai's model as 
well equilibrated well which get emptied on aging time scales. This picture should be checked further. 

Another consequence of our result for Sinai's model is that, since the aging regime must be compatible with the 
diffusion one and we have found that h{t) ~ L(t) ~ Int, there must be a singularity in the two time diffusion front at 
2: = to allow for a non trivial aging regime. 

We have given a general explanation of these regimes using scaling arguments on the next highest barrier encountered 
by the particle. This allows to understand the aging form in Sinai's model. It also strongly suggests that the aging 
in Sinai's model could be studied analytically by only computing the distribution of next highest barrier (a purely 
geometrical feature of the energy landscape) . These considerations also lead to define a class of models for which the 
distribution of these barriers can be computed exactly, and allows for predictions of the aging forms (159). 

Though we did get a consistent picture of aging in Sinai's model, we cannot rule out completely other regimes. For 
instance, we have not explored in more details the behaviour of the moments of the displacement in Sinai's model. 
As in Section III one could say that at the very beginning of the aging regime (i.e Int/lnt^ ^ 1 + e fixed, a small 
fraction of particles have escaped from their well and have experienced Sinai's diffusion to another well. One then 
gets: 

I < x{t) - x{t^) > I" ~ (-^ - lf{\n{t - i^))2" ~ (lnt^)2"-/3(lnt - Int^f (161) 
In t-iu 

Thus by the same mechanism as in Section III the various moments may have some different aging behaviours. We 
have not attempted to obtain a precise estimate for /3 but a rough estimate from our numerical simulations (Fig. 13 
) is consistent with (3 = 1. If this is the case the moment n = 1/2 may have an aging behaviour as t/tw The general 
issue of the matching bewteen the three regimes defined here deserves to be investigated further. 

Another open problem is the behaviour of the response in these regimes. It is important to determine how to define 
properly, and study beyond mean field the way the equilibrium theorems are violated. Since, as we have shown, 
sample to sample fluctuations play a strong role and one should focus on distributions, one needs extensions of the 
mean fleld ideas. As a flrst step we have given analytical expressions and deflnitions of quantities adapted to low 
dimension and which measure these violations. A detailed numerical and further analytical investigation of these 
quantities is deferred to the future. 

To summarize, we have found that some of the concepts deflned in mean fleld are still useful in low dimensional 
models, though they have to be seriously adapted. We hope that this study will also help understand dynamical 
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behaviour in low dimensional but more complex systems such as domain wall motion with disorder and coarsening in 
random spin systems, where ultra slow anomalous diffusion processes are expected to play a crucial role. We expect 
that the various regimes defined here should be present in these systems as well. 

We thank A. Barrat, A. Georges, L. Cugliandolo and J. Kurchan for useful discussions. 



APPENDIX A: DISCRETE VERSION OF SINAI MODEL AND DIAGONALIZATION 

In this Appendix we describe in detail the observables for discrete hopping models. We also describe the numerical 
method used in the paper. We establish some FDT relations and other useful exact relations for discrete models. 
Let us consider the Fokker Planck operator Hpp defined in (7). It can be written as: 

dPn 

{HFp)n,mPm = = —{Jn+l,n — Jn,n-l) (Al) 

where the current flowing from site n — 1 to n is by definition: 

Jn^n-i = - e-^^P™ (A2) 

We need to compute the Green function P{n,t\no,to) {t > to) which is defined as the solution of (7) with initial 
condition P{n,to\no,to) = Snno- Following Ref.^® it is useful to map the FP equation onto a Schrodinger equation 
corresponding to a symmetric matrix. One has: 

P{n,t\no,to) = e-iC^-^-o) ^ v^^^^oe""""* (A3) 

a 

where the are the eigenstates of the Schrodinger operator: 

= -(C+l + C-l - 2C) + = Eai^n (A4) 

in the potential: 

Vn = e'^-'+i + e-^" - 2 (A5) 

As discussed in Ref.^^ the random operator Hg, which is a version of supersymmetric quantum mechanics, is quite 
peculiar : all states are localized but the spectrum of Hs is positive, without the Lifschitz tails usually associated 

to random one-dimensional potentials. An eigenfunction corresponding to the energy level ii'o = is always exactly 
known, i.e ^/;° — Ze^^'^l'^ . Whether or not this is the actual ground state depends on whether -^J^ is normalizable, i.e 
on the boundary conditions. This is also related to the breaking of supersymmetry. 

The Schrodinger operator is a tridiagonal symmetric matrix and is easily diagonalized for large size L. Let \is first 
consider the problem with L + 1 sites k = 0,L with reflexive boundaries. The same change of function P„ — e^^/'^tj^n 
is used and formula (A3) holds. 

{HFp)oj = Sije-'^' - Soje'^' 

{H,)oj = -Sij + Soje*' (A6) 

these boundary conditions simply amount to choose 4>-i = +oo and 4>l+i = — oo and restrict the problem to sites 
k = 0,..L. 

explicit expressions of quantities of interest and FDT theorems in discrete version 
The two time quantities of interest are correlation functions of some operator 0{n,n'): 

0{t,t') = < 0{x{t),x{t')) >=^0{n,n')P{n,t\n',t')P{n',t'\nQ,tQ = 0) 

n.n' 

= E e-^<^(*-*')-^^*' E 0{n, n')e-^"/^CC'^f 'e^-/2^^^ (A7) 

a,/3 n n' 
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and response functions of some operator R{n, n') 

i?(M') = (A8) 

defined by adding a short-duration pulse of small uniform force of integrated strength /, i.e <pn (j>n + f at time 
t'. Using the definition P{n,t\n',t') = (e^*'^'(*"*'))„„' one has: 

R^{t,t') = l Yl R{n,n')P{n,t\n',t')Br,,rn'P{rn',t'\no,to = Q) (A9) 

if the pulse is at t' — e, or 

RM = \ R{n,n')P{n,t\m,t')Bmn'P{n',t'\no,to = 0) (AlO) 



2 

if the pulse is at t' + e (Ito). 
We have defined 

J^mm' — / J ri — <Jm,Tn'('i e ) ^ tJm—l,m''i Om-|-l,m'e 

0(Pn' 

n 

We will now establish (i) the FDT relations valid when equilibrium is attained (ii) some exact relations always valid. 
These equations will relate the correlation functions of the operator 0{n,n') with response function of the operator 
R{n,n'). There must be a relation between R and O for these to hold, which is: 

i(i?(n, m) + R{n, m - 1)) = 0(n, m) - 0{n, m - 1) (A12) 

This relation generalizes the usual relation R{x,x') = dx'0{x,x') (see next Section) valid for continuous systems. 
We start with the following identity: 

^£M1= 0{n,n')P{n,t\m,t')Hmm'{Smn' -Sm'n')P{m',t'\no,to = 0) (A13) 

Note the simplification: 

Hmm' {.^mn' ^m'n'^ ~ i^^m+l,m'^ ~t~ ^m—l^m'^^^ '^^^{.^mn' ^m'n') (A14) 

Now the following exact relation can be established: 

0{n,n')Hmm'{6mn' - 5m'n')Pm' = ^ nQ-Bw f +\R{n,m){Jm+l,m + Jm,m-l) (A15) 

n' ,m' n' 

for any set of Pm, provided the above relation (A12) holds between the operators O and R (J being defined as in 
(A2)). It yields to: 

^^gf^ = R{t,t') + i^i?(n,m)P(n,t|m,t')(^n^+i,m(i') + J^,m-i{t')) (A16) 

where Jm,m-i{t') = e'^"'P{m — 1, t'\no, to = 0) — e~'^^P{m, t'|no, io = 0) and we are using the Ito response. 

When equilibrium is attained, i.e either in the limit t' +oo before L +oo or, if there is an FDT regime in the 
problem (see text, this usually entails averaging all these correlations over disorder), then one can set the current to 
zero, J = 0. Then the the fluctuation-dissipation relation holds: 

^ = K(*,0 (A17) 

applications 
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- The choice 0{m,n) = nm and R{n,m) = n is consistent with (A12). It gives: 

d < xit)x{t') > 6 < x{t) >f 

— dt' — = sm ^^^^> 

- The choice 0(m, n) = (n — m)^ and R{n, m) = — 2(n — m) is consistent with (A12). It gives: 

d<{x{t)-x{t')r> _ j<{x{t)-x{t') >f 
dt' ~ ^ sm 

- The general choice: 

0{n,m) = ^{d{n,m)+d{n,m+l)) R{n,m) = d{n,m + 1) - d{n,m) (A20) 
satisfies the condition (A12). 

Thus one can choose 0{n,m) = and obtain the exact relation: 

dt'^iQiz, t, t') + Q{z + 1, t, t')) = R{z + l,t, t') - R{z, t, t') (A21) 
+ 1 ^(^('^ + 1 + m, t\mt') - P{z + m, t\mt')){Jm+i,m{t') + Jm,m-i{t')) (A22) 

m 

This is the discrete equivalent of the continuous relation derived in the next Section (B23). 

Finally, one can wonder what happens when a bias is applied on a finite size periodic ring. There a stationary 
distribution with a fixed current is reached at large t' . There is an extension of the FDT theorem. Indeed one has: 

^^1^ = R{t, t')+jY^ R{n, m)P{n, t\m, t') (A23) 
Note that J is simply related to the velocity (see e.g^^). 

APPENDIX B: FDT AND USEFUL EXACT RELATIONS FOR PROBABILITY DISTRIBUTIONS 

In this Appendix we derive a generalization of the FDT theorem on the probability distribution. It can then be 
used as a generating functional to obtain a hierarchy of FDT relations on all moments of the type < a;(i)"a;(t')'" >. 
We are interested in the joint probability that the particle is in x at t' and then in x at t. 

P{xt, x't'\xoto) = P{x, t\x',t')P{x', t'\xo,0) (Bl) 

Let us recall the forward and backward FP equations: 

dtPixt\x't') = Td^D(x)d^P{xt\x't') - d^D{x)F{x)P{xt\x't') (B2) 
dfP{xt\x't') = -Td^>D{x')d^'P{xt\x't') - D{x')F{x')d^.P{xt\x't') (B3) 

We now derive the differential equation for the joint probability P: 

dt'P = {dt'P{x, t\x', t'))P{x', t'% 0) + P{x, t\x', t'){dt'P{x', t'\0, 0)) (B4) 

dt'P = -J dydy'P{x, t\y, t'){[HFP, 6,,])yy,P{y' , t'\0, 0) (B5) 

with: 

[Hpp, S,,] = [TdDd - dDF, S,,] = TdD[d, 6,,] + [8, 5,,]{TDd - DF) (B6) 

The last term is the current J = —{TDd — DF)P. In the FDT regime (for large t') the current is expected to 
vanish and we are left with the following equation for P: 
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9t'P = -J dydy'Pix, t\y, t'){TdD[d, S,,])yy,Piy' , t'\0, 0) (B7) 

dt'P = d^, J dydy'P{x,t\y,t'){TdD5,,,)yy,P{y',t'\0,0) (B8) 

dt>P = -Td^.[{d^,P{x,t\x',t'))D{x')P{x',t'\0,0)] (B9) 

finally we obtain: 

dt'P = -Tdl,{D{x')P) +Td^,{P{x,t\x' ,t')d^,{D{x')P{x' ,t'\Q,Q))) (BIO) 
and since = —dx'D{x') 

dt>P = -Tdl,{D{x')P) - (Bll) 

where we define P^{xt\x't'\xoQ) the joint probability when a field pulse has been applied at t' — e. This is the 
equation which relates exactly the joint probability distribution P{xt,x't'\xf)tQ) to the response distribution in the 
quasi-cquilibriuni FDT regime. It was obtained by setting the current at time t' to zero. 

A similar equation can be derived for the field applied at time t' + e. From (BIO) one has also that: 

^*'^ = -^^-^ ^^''^ 

where Ph{xt\x't' \xqO) the joint probability when a field pulse has been applied at t' + e. This corresponds to Ito's 
prescription for the response functions since SPh{xt'\x't'\xoO)/6h{t') = 0. 

Prom this it is immediate to derive a similar FDT equation for the probability function Q{z,t,t') = J dxdx'6{z — 
X + x')P{xt,x't'\xf)tf)). Multiplying the above equation by the delta function, integrating with respect to x and x', 
and integrating by parts one gets: 

dt'Qiz, t, t') = -TdlQ{z, t, t') + Td, ^^I'f'/j; + Td, [ dxidx2S{z-xi+X2)P{xi,t\x2,t')J{x2) (B13) 

dh{t') J 

in the the long time regime — > oo the current vanishes and we are left with the FDT regime of the equation: 

dt,Q{z,t,t') = -TdlQ{z,t,t')+Tdf-^^^0^ (B14) 



Similarly with the Ito prescription: 



d,Q{z,t,t')=Td.^-^^^^^ (B15) 



We can check that this equation is more general than the conventional FDT theorem, and indeed gives back the 
usual FDT result. Defining: 

B(t,t')FDT = j dzz^QFDT{z,t,t') (B16) 
and inserting it in the generalized equation gives: 

dt^B{t,t') = -T j dzz^dlQ{z,t,t')+T j dzz'dj-^^^^ (BIT) 
after integration by part this gives simply: 

dt'B{t,t') = -2T-2Tjj^ J dzzQ+{z,t,t') = -2T{1 + R+{t,t') - R+{t',t' - e)) = -2TR{t,t') (B18) 
With the Ito prescription one has simply: 
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dt'B{t, t') = -2T^^ j dzzQhiz, t, t') = -2TR{t, t') (B19) 
using R{t', t') = R{t', t' + e)= 0. 

A motivation is to find a generalized form to this equation which would be valid in the aging regime as well. 
exact relations and FDT violation ratios 

It is useful also to give the exact relations (always valid) for averages of operators. They allow to obtain explicit 
the FDT violation ratios. 

We put back the current term that we have neglected. We obtain then instead of (BIO) the (still) exact relation: 

dt'P + Tdl,{D{x')P) - Td^,{P{x,t\y,t')d^,D{x')P{x' ,t'\Q,Q)) = j dydy' P{x,t\y,t')[d,5^,]yy,J{y't'\m) (B20) 
with J{y't'\m) = -D{y'){Tdy' - F{y' j)P{y't'\00). We wiU use the Ito response here. 

Let us study a general observable 0{x, x'). Multiplying the above equation, integrating over x and x', and integrating 
by parts (assuming no contributions from boundaries) one has an exact relation which relates the averages of of 0(a;, x') 
and of 0'{x,x') = dx>0{x,x'). Defining: 

Vo{t,t') = dt, < 0{x{t),x{t')) > -T^^^^l^^ = (1 - Xo{t,t'))dt' <0> {t,t') (B21) 



the relation reads: 



Vo{t,t')= j dxdx'P{x,t\x',t')0\x,x')J{x't'\m) (B22) 

We have defined above the FDT violation ratio Xq associated with the operator O. If the other response was used 
there would be in addition a term — T < D{x{t'))02{x{t),x{t')) > in the above equation with 02{x, x') = d^,0{x, x'). 
Again for the Ito response the second derivative term is absent. 

We will also give an exact relation for the relative displacements. Let us consider an operator 0{z). Then one has: 

dt, < 0{Z) > +r '^^£}jj = fdXidX20'{Xi - X2)P{Xit\X2t')J{X2t'\00) (B23) 

oh{t') J 
And thus one can define also a generalized FDT ratio: 

~ JdxidX20'{xi - X2)Pixit\x2t')Jix2t'm 

^oit,t) = i a,, < o{z) > (^^^^ 



APPENDIX C: BOUNDS 

Here wc illustrate the bounds recently proposed by CDK'*" . We use the framework of the generalized FDT relation 
of the preceding section and our derivation is thus technically slightly different, though identical in spirit to"*". We 
work directly with the FP equation and a space dependent diffusion coefficient. 

The nice observation of CDK is that the current which appear in (B22) also appears in the H theorem which states 
that the free energy: 

H{t') = J dx'P{x't'\00){TlnP{x',t'\00) - U{x')) (CI) 

is always decreasing with: 

dim _ _ f , , (j(.x't'ioo))^ 

dt' ~ j D{x')p{x',t'm ^ ' 

The CKS bounding amounts to bound Vo{t,t') defined in (B21) by: 
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\Vo{t,t')\ < l^^^r^^y dxdx'P{x,t\x',t')O'{x,x'fD{x')P{x',t'\00) (C3) 
= \^\'/' < {0'{x{t),x{t')fD{x{t')) (C4) 
using the Cauchy-Schwartz CS inequality and J dxP{x,t\x't') — 1 

Note also that disorder averages can be bounded similarly by applying the CS inequality at the same time to the 
integrals over x, x' and configurations. 
One then gets: 

mm\ < i^^r/'i< o'{xit),xit')fDix{t')) >\'/' (C5) 

Or, as pointd out by CDK, in an integrated version: 

\j\sV^)\ < £ ds{\^^\'^^\< 0'{x{t),x{s)rD{x{s)) >\'/^) (C6) 

If we choose 0{x, x') = xx' and D{x) = 1 one gets: 

\< x^t) > - < x{t)x{t') > - T f R{t,t')\ < (< x^t) >)i/2 /*rf5(|^^|i/2 (C7) 

Jt' Jf as 

One can also derive bounds using (B23) for the relative displacements. Using CS, (B23) leads to the bound: 

\d, < 0{z) > +T^^^[||^| < I < {0'{x{t) - x{t')f > r/^l^r/^ (C8) 
Which can be rewritten as: 

\d, J dzO{z)Q{z,t,t')+T^^ J dzO\z)^-^^^01\ < I J dzO'Hz)Q{z,t,t')n^\'/' (C9) 
This yields in particular: 

\dt'B{t,t') + 2TR{t,t')\ < B(t,t')V2|^^|i/2 (CIO) 

Or its integrated version: 

\-B{t,t')+2T f R{t,t')\< f dsB{t,sf'^\'^^^\^'^ (Cll) 
Jt' Jt' ds 

More generally the bound can be used to constrain the Xq defined in (B24): 

< 0'2 > (t 



APPENDIX D: MAPPINGS OF SEVERAL MODELS 

The method of change of variables allows to relate exactly members of a class of landscape. Let two landscapes 
and their corresponding Green's functions be: 

{Eb{x), U{x)) ^ P{x, t\xQ, 0) (Dl) 
{E'i,{x),U'{x))^P'{x,t\xo,0) (D2) 



If there exist a function y{x) such that: 
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E',{x)=E,{y{x)) + 2ln^ U'{x) = U{y{x)) -In^ (D3) 
then the two Green functions are related through: 

P'{x, t\xo, 0) = ^P{y{x),t\y{xo), 0) (D4) 

In particular one can map: 

(0,C/(a;))^(2[/(x(w)),0) (D5) 

and 

{Et,{x),0) ^ {0,^Eb{x{u))) (D6) 

It is interesting in general because the new functions U{x{u)) or Eb{x{u)) are usually better behaved at large u. 

Two general scenarios exist, confining potentials, unconfining ones. Let us take rapidly growing landscapes {Eb{x), 
U{x)) as in last section and map them onto (0, U'{u)). Then typically one has U'{u) ~ h\nu at large u. The case 60 
is confining and corresponds to the case where barriers grow faster than valleys Eb{x) — 2U{x). The case 60 is fast 
diffusion and corresponds to the case where barriers grow slower than valleys Eb{x) < —2U{x). 

Finally note that time dependent mappings could also be built using the following propagator: 



which satisfies: 



dtQ{u, t) = duduQ{u, t) - ^duuQ{u, t) (D8) 



APPENDIX E: RESPONSE IN THE AGING MODEL 



We start from the exact equation obeyed by the Brownian diffusion propagator: 

1 

R{ut\u't') = =e ^-^c-*') (El) 

V47rr(t - t') 

The joint propagator 

R(ut\u't'\uotQ) = R{ut\u't')R{u't'\uoto) (E2) 

satisfies the exact equation: 

dt'R = -Tdl.R + 2Tdu-{R{ut\u't')du'R{u't'\uotQ)) (E3) 

This is valid for the free (unbounded) brownian motion and yields for instance R(t,t') — Xdt'C{t,t') with X = 
1/(2T). A bounded brownian motion would instead converge to equilibrium with X = 1/T and satisfy the FDT 
equation with 2T replaced by T in the last term. 

One can now use the change of variable ^{x)d/dx = d/du and R{ut\u't'\uoto) = ^{x)^{x')P{xt\x't'\xoto) and 
obtain for the model studied previously, the exact equation valid in the unbounded case: 

dt,P = -Tdl,{<^{x'YP) + Td,.{P{xt\x't')d,.{^{x'fP{x't'\xoto))) (E5) 
+Td^,{P{xt\x't')^{x')d^,{^{x')P{x't'\xoto))) (E6) 

Since the diffusion coefficient is D{x') — (f>(a;')^ in this model, and the response is |^f7j = —dx'D{x') the above 
equation can be rewritten in a form very similar - but not identical - to the above general FDT equation: 
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SP, 
5h{f') 



'--d^,{P{xt\x't')J{x't')) 



(E7) 



The last term cannot be simplified further and involves the current J{x't') = —T(^{x')dx'{f^{x')P{x't'\xQtQ)) of the 
model. This equation is always valid for our model, even in the out of equilibrium regime. 

The idea is that in the non trivial aging regime all three terms of the above equation will be roughly of the same 
order in t' and thus it will effectively lead to a non trivial FDT ratio X{t, t'). 

One can use e.g this equation to study the correlation C{t,t') =< x{t)x{t') >. Multiplying by xx' and integrating 
over X and x' one obtains: 



= ^ ^Jin = X{t,t')dt' < x{t)x{t') > 



where we have defined: 



One can also define: 



Sh{t') 



X{t,t') = -{1- 



J dxdx' xP{xt\x't')J{x't') 
df < x{t)x{t') > 

dt,<x{t'f 



Introducing: 



A{t',T) =< xlVVui + ^U2]ui<^[x[\/t'Ui]] 

D{t',T) =< (iti - U2^)x[Vt'Ui]^x[Vt'Ui 



'TU2\ 



One finds: 



Xc = 

Xb = 



D 



TA + D 
1 D 



2TD + A-A0 



where Aq = A{t',0). 

Thus one gets for instance: 



(E8) 
(E9) 

(ElO) 

(Ell) 
(E12) 

(E13) 

(E14) 



Xb 



X < e(ui)e(ui + XU2) ln(l + Vt'lm + xu2\)~ < ln(l + Vt'\ui\) 



2T 



< f7T^e("i)<"i + 2^"2) ln(l + Vf\ui + xu2\) 



(E15) 



with X = . The divergences make the calculation depend strongly on the boundary conditions. As explained 
in the text we have not pursued it further. 



Defining: 



We will use that: 



APPENDIX F: DIRECTED MODEL CALCULATIONS 



W 



1 - s$(s) 



1 $(.Si)-$(s2) 



(Sl + W){S2 + W) S2 - Si 

W _ S2HS2) - si$(si) 



{SI+W){S2 + W) S2-S1 
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Let us compute the averaged probability that the particle advances by m between t' and t = t' + t: 
Q{m, T, t') = {S{x{t) - x{t') - m)) = ^ P{n + m,n,T)P{n,Q,t') 



n>0 



The double LT, Q(m,si,S2) = drdt'e-^^''-^^*' P{m,T,t') can be calculated: 



Q{m, si, S2) = V P{n + m, n, si)P{n, 0, S2) = {Smo{j -^77^7 -^777) 

^ {S1+W){S2 + W) 



It yields the result given in the text. 

Let us now estimate the probability Pt^ (W) that at time the walker is on a site with a waiting time W 
Its Laplace transform with respect to t,,, (Laplace variable .S2) is simply given by: 



n k=no 



This easily leads to: 



S2 ^S2 + W^{S2) 

By Laplace inversion this yields for the distribution of the waiting time f = 1/W: 



TT tw T Jq T[ij,\ 
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